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ABSTRACT 


This  volume  summarizes  that  portion  of  the  program  concerned  with  the 
development  of  theoretical  methods  of*  analysis  for  orthotropic  laminated 
conposite  structure.  An  appendix  to  this  volume  presents  a  nunber  of  conputer 
programs  that  were  generated  in  conjunction  with  tiiese  theoretical  developments. 

The  main  bocfy  of  this  volume  is  comprised  of  independent  sections  devoted 
to  the  following  topics: 

Extensional  and  flexural  elastic  constants 

Flat  plates 

Orthotropic  .'sandwich  plates  and  shells 

Stiffened  skin  construction 

Minimum  wei^t  analysis 

Thermoelastic  relationships 

Joint  analysis 

The  appendix  to  this  volume  presents,  for  each  of  the  included  conputer 
programs,  an  abstract,  data  input  instructions,  output  description,  and  a 
typical  sxanple  problem.  The  listed  programs  cover  the  following  subjects: 

Calculation  of  laminate  extensional  and  flexural  elastic  constants. 

Honeycomb  sandwich  panels  under  in-plane  biaxial  loading,  in-plane 
shear  loading,  nonnal  pressure,  or  certain  combinations  of  these. 

Hat-  and  Z-stiffened  panel  stability  and  stiffness. 

Laminate  thermoelastic  analysis. 

This  abstract  is  subject  to  special  e:qx)rt  controls, 
and  each  transmittal  to  foreign  governments  or 
foreign  nationals  may  be  made  only  with  prior  ^proval 
of  the  Air  Force  Materials  Laboratory  (AIM./LC) , 

Wright-Patterson  Air  Force  Base,  Ohio  45433. 
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Section  I 


INTRODUCTION 


The  purpose  o£  this  program  was  to  take  the  first  step  toward  the 
generation  and  presentation  of  basic  engineering  data  necessary  to  perform 
high-confidence-level  structural  design  of  primary  aircraft  structures  utilizing 
advanced  composite  materials.  The  program  was  limited  to  an  in-depth  generation 
of  basic  material  allowables  for  one  boron/epoxy  and  several  graphite/ epoxy 
material  systems,  and  the  determination  of  basic  structural  element  response 
for  the  boron/epoxy  system  alone.  The  boron  portion  of  this  program  was  con¬ 
ducted  in  conjunction  with  a  concurrent  General  Dynamics/Forth  Worth  (GD/FIV) 
program  which  was  funded  under  Air  Force  Contract  F33615-68-C-1474.  Tlie 
boron/epoxy  material  system  highlighted  by  both  these  programs  was  Narmco  5505, 
furnished  by  the  supplier  as  3-inch  prepreg  tape.  This  is  a  composite  material 
consisting  of  collimated  4 -mil  boron  filaments,  208  per  incli  of  tape  widtli, 
embedded  in  a  matrix  of  Narmco  2387  epoxy  resin,  and  supported  on  a  1-mil  layer 
of  104  glass  scrim  cloth.  The  graphite  portion  of  this  program  is  being  con¬ 
ducted  independent  of  any  other  program,  I  will  be  described  in  fuller  detail 
in  a  later  volume  of  this  report.  AdditiOiiul  data  for  all  these  materials,  as 
well  as  for  other  filament/matrix  material  systems,  were  obtained  from  published 
Government,  industry,  and  technical  journal  reports,  and  were  used  to  augment 
the  data  generated  in  this  program. 

This  program  was  composed  of  three  major  work  task  areas: 

Task  I  -  Generation  of  Composite  Material  Design  Allowables 

Task  II  -  Structural  Element  Test  Program  and  Analysis  Evaluation 

Task  III  -  Development  of  Advanced  Con^josite  Structural  Design  Manual 
for  Aircraft 

Task  I  is  divided  into  two  distinct  areas  of  effort  by  the  separate  boron/ 
epoxy  and  graphite/epoxy  programs.  The  purpose  of  the  boron  portion  of  task  I 
was  to  complement  the  basic  material  design  allowable  activities  conducted  by 
GD/FIV  (reference  12)  and  to  develop  acceptable  laminate  fabrication  and  inspec¬ 
tion  procedures.  The  boron  effort  was  divided  into  the  following  work  areas: 
Tlie  establishment  of  program  coordination  procedures  for  the  North  American 
Rockwell  Corporation  and  General  Dynamics  related  programs;  the  acconplishment 
of  a  limited  material  development  program;  the  generation  of  basic  allowables 
for  the  constituent  materials;  establishment  of  the  accuracy  of  current  analy¬ 
tical  procedures  for  predicting  certain  basic  allowables;  and  the  development, 
where  reliable  techniques  were  lacking,  of  prediction  techniques  for  these 
basic  material  allowables. 
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Hie  graphite  portion  of  task  I  consists  primarily  of  a  screening  and 
characterization  of  several  graphite/epoxy  material  systems  nnd  will  be  delin¬ 
eated  in  a  later  volume  of  this  report. 

Hie  purpose  of  task  II,  which  was  concerned  solely  with  boron/epoxy 
material,  was  to  generate  data  on  basic  structural  elements  whicli  form  tiie 
building  blocks  from  which  aircraft  structures  are  designed.  A  minimum  evalua¬ 
tion  of  structural  elements  was  conducted,  including  one  basic  laminate  and 
one  elevated  temperature.  Factors  which  were  considered  in  the  detail  design 
of  the  structural  elements  included  laminate  orientations,  panel  proportions 
and  edge  restraints,  effectiveness  of  typical  forms  of  panel  stabilization, 
evaluation  of  cutouts,  and  thermal  gradient  effects.  One  or  more  elements 
were  selected  for  each  primary  and/or  combined  load  applications.  Hie  test 
program  included  local  and  general  instability  of  flat  panels  and  natural 
frequency  determinations.  Hie  results  of  this  test  program  were  compared  to 
predicted  response,  failure  mode,  and  strength  techniques  for  basic  structural 
elements . 

The  task  III  work  area  was  originally  centered  on  the  development  of  an 
advanced  conposite  structural  design  manual  for  aircraft  structures.  The  first 
effort  of  this  task  involved  revision  and  refinement  of  the  Aircraft  Division 
of  the  Intermediate  Draft  of  the  Structural  Design  Guide  developed  by  the 
Southwest  Research  Institute,  San  Antonio,  Texas,  under  Air  Force  Contract 
AF33  (615) -5142.  The  completely  revised  and  reorganized  Aircraft  Division 
resulting  from  this  phase  of  effort  was  published  in  the  Final  Draft  of  the 
Design  Guide  in  November  1968  under  Contract  AF33(615)-68-C-1241.  Soon  there¬ 
after,  a  review  of  the  Final  Draft  by  a  select  industry  group  led  to  a  decision 
by  AFML  to  reorganize  the  entire  Design  Guide  for  the  First  Edition,  which  was 
then  assigned  to  NR/LAD  under  Contract  F33615-69-C-1368. 

Subsequent  phases  of  task  III  of  this  program,  in  light  of  the  foregoing 
developments,  consisted  of  the  preparation  of  the  Aircraft  System  Applications 
chapter  of  the  First  Edition  of  the  Design  Guide  as  well  as  the  preparation  of 
data  generated  by  tasks  I  and  II  of  this  program  for  incorporation  into  the 
various  technical  function-oriented  chapters.  Task  III  also  included  the 
incorporation  into  the  Design  Guide  of  data  generated  by  the  concurrent  GD 
program. 

The  bulk  of  the  basic  material  allowables  for  the  5505  material  system  was 
generated  by  the  General  Dynamics  contract.  This  concurrent  and  integrally 
related  contract  was  coordinated  with  the  Los  Angeles  Division  program  effort 
through  scheduled  periodic  coordination  meetings.  These  meetings  insured  the 
continuous  flow  of  pertinent  program  data  between  the  two  contractors. 

This  report  is  divided  into  four  separate  volumes,  in  each  of  which  the 
subject  areas  of  interest  comprise  an  independent  segment  of  the  overall  pro¬ 
gram.  Each  volume  is  a  self-contained  document,  conplementing  the  other  three 
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volumes  but  not  dependent  upon  them  for  coherence  or  continuity.  The  titles  of 
the  four  volumes  are: 

Volume  I  -  Material  and  Basic  Allowable  Development  -  Boron/Epoxy 

Volume  II  -  Structural  Element  Behavior  -  Test  and  Analytical 
Determination 

Volume  III  -  Tlieoretical  Methods 

Volume  IV  -  Material  and  Basic  Allowable  Development  -  Graphite/Epoxy 

Volume  II  contains  the  analytical  derivation  of  the  behavioral  equations 
and  relationships  characteristic  of  anisotropic  composites.  It  is  divided  into 
a  number  of  independent  sections,  each  of  which  treats  a  particular  structural 
element  concept  or  some  separate  facet  of  advanced  composites  application  to 
structural  design. 

Section  II  is  fundamental  in  scope  and  concerns  the  extensional  and  flexu¬ 
ral  elastic  coiistants  of  filamentary  composite  laminae  and  laminates.  Sections 
III,  IV,  and  V  cover  the  behavioral  characteristics  of  flat  plates,  orthotropic 
sandwich  plates  and  shells,  and  stiffened  skin  construction,  respectively. 
Section  IV  discusses  the  general  approach  to  minimum  weight  design,  with 
detailed  attention  to  zee-  and  hat-stiffened  elements,  while  Section  VII  is 
devoted  to  the  thermoelastic  relationships  of  composites.  Section  VIII  con¬ 
cerns  the  special  topic  of  con^^osite  panel  joints,  and  includes  separate 
treatment  of  several  types  of  joint. 

Laminate  ply  orientations  are  described  and  specified  in  this  report  by  use 
of  the  laminate  orientation  code  defined  in  the  Structural  Design  Guide  for 
Advanced  Composite  Applications. 
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SECTION  II 


EXTENSIONS  AND  FLEXURAL  ELASTIC  CONSTANTS 

STRESS-STRAIN  RELATIONS  OF  A  LAMINA  IN  GENERALIZED  PLANE  STRESS 


Consider  the  fiber  reinforcement  to  be  unidirectional.  Referring  to  a 
set  of  orthogonal  axes,  L,  T,  and  Z,  with  L  parallel  to  the  fiber  direction 
and  Z  =  0  coinciding  with  the  middle  surface  of  the  lamina,  then  the  stress - 
strain  relation  of  the  lamina  shown  in  figure  1  may  be  represented  by 
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where  <tl»  ^T  direct  stresses,  ffi;i  is  the  shear  stress,  ei,  ej  are  direct 
strains,  andy^T  shear  strain. 
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Figure  1.  Typical  Lamina 
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El>  Et  denote  the  Young's  modulus  in  L  and  T  directions,  respectively, 
denotes  the  shear  nwdulus ,  and  I'XL  denote  the  major  and  minor 
Poisson's  ratios,  respectively. 

Let  X,  Y,  and  Z  be  a  set  of  rectangjilar  axes  vdiere  the  X  and  Y  axes  are 
obtained  by  rotating  L  and  T  axes  throu^  an  angle  6  as  shown  in  figure  2. 


Figure  2.  Lamina  and  Laminate  Axes 

The  stress-strain  relation  referring  to  X,  Y  axes  is  obtained  by  transforming 
equation  1  as 


11  ^12  ^161  ^X 


^Y  '  ^21  ^22  ^26  '  ^Y 
‘'XY  ^61  ^2  ^66 


where 


Qn  =  <511  •  «22 


*^12  ~  *  '^22  '  ^*^6^  ”  ”  *  ‘^12^’"  *  "  ^  "  '^21 


•Jw  ■  «n  -  ‘^12  -  ”  "  *  «12  -  ^22  *  2Qee)  *  = 

^22-'  “’ll  ""  2^12  ^  *  '^22 
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‘^6  '  «11  "  '^22  -  2Qi2  -  ^6’  *  %6 


m  =  cosO 


n  =  sinO 
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It  should  be  noted  that  the  stress-strain  relation  referring  to  the 
laniina  axes  (equation  1)  is  orthotropic,  i.e.»  no  coupling  between  the  exten- 
sional  and  shearing  stresses  and  strains,  while  the  stress-strain  relation 
referring  to  the  laminate  axes  (equation  2)  is  generally  not  orthotropic. 


SECTIONAL  FORCES  AND  MOMENTS  OF  A  LAMINATE  IN  TERMS  OF  DlSPLACBenS 

From  Kirchhoff*s  assunption,  normals  to  the  middle  plane  of  the  laminate 
remain  normal  to  the  middle  surface  of  the  l.miinate  during  deformation,  and 
the  normal  stress  on  planes  parallel  to  the  middle  surface  may  be  neglected 
in  conparison  with  other  stress  conponents.  Let  u,  v  and  w  be  the 
displacement  along  the  X,  Y  and  Z  axes,  and  let  Uq,  vq  and  vjq  be  the 
corresponding  displacement  of  points  on  the  middle  plane.  We  have 


u  =  u-  -  z  s~ 
0  dx 

Sw 

V  =  V-  -  z  5- 
0  dy 


w  =  w^ 


For  small  strain. 


X  dx*  y  dy*  ^xy  \dy  dx) 
Substituting  equation  4  into  equation  2,  we  have 


ryJi  sy  a/ 
^ex  ay  / 


d30y\ 


where  the  subscript  "i”  denotes  the  i^  lamina  in  the  laminate,  TTie  notations 
of  the  laminate  are  shown  in  figure  3.  For  a  laminate  of  N  laminae,  the 
section  forces  and  mcments  shewn  in  figure  4  are  defined  as 
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Figure  3.  Laminate  of  N  Laminae 


Fig;uTe  4.  Positive  Direction  of  Internal  Forces  and  Moments 


The  sign  conventicns  for  those  section  forces  and  moments  are  shown  in  figure 
4.  The  laminae  are  now  considered  to  be  synmetrically  laid  vdth  respect  to 
the  middle  plane  of  the  laminate.  Perfoiming  the  integration  as  shewn  in 
equation  6,  we  have 


Let  us  now  further  consider  only  the  cases  where  the  number  of  laminae 
with  +6  is  equal  to  that  with  -0,  i.e. ,  the  layup  is  balanced.  Then 
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and  equation  7a  reduces  to 
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The  stress-strain  relation  corresponding  to  equation  11  for  a  homogeneous 
plate  is  given  as 

S  ^11  ^12  ° 

'  <^y  ' "  t  \l  ^22  °  '  V  ' 
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Comparing  this  with  equation  1,  it  is  seen  that  the  equivalent  homogeneous 
elastic  constants  of  the  laminate  are  then 

V  ' 

V  '  V''22 


^1  ^22  ■  ^12 


A  A  -  A 
11  22  12 


xy  t 


In  addition,  since  the  laminate  is  considered  to  be  laid  up  symmetrically 
relative  to  its  midsurface,  equation  9  reduces  to 


D^^'s  OF  SOME  COMMONLY  USED  LAMINATES 

Let  us  now  consider  symmetric  laminates  somposed  of  q  sets  of  plies  with 
each  set  having  p^  plies  at  0°,  2p^  plies  at  ±0° ,  and  po  plies  at  90°  as  shoi-m 
in  figure  5. 

Let 


Figure  5.  Typical  Synmetrical  Laminate 
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Letting 


2  ■  ^ 
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^3  =  2  -  P2 
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Expanding  <p  in  equation  18  and  noting  that 


n^  =  I  (N+l)  (2N+1) 
n=l 


E  "  =  f  (N*i) 
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Cj^  "  -4-^  CPj+1)  [2Pj  +  H-  3  (2rj  -  N  -  1)1  +  |  (2rj  -  N  -  1)  +  j 


(3-*  1) 

+  3ipj  (pj  »  2rj  -N)  a 


^*3p  P^adlilH 
2  2  6 


q  Pi  r  2 

=  -^[2(p^+l)  C2Pi  +  6r^  ■■  3N  -  2)  +  3  C2r^  -  N  -  1) 

+  r  (q+2)  (3p^  +  6r^  -  3N  +  r  +  qr)  +  ij 

By  replacing  by  and  p^^  by  2p^  in  the  preceding  equation,  we  obtain 

r  2 

S  "  ~T  |_2(2p/l)  (4p^+6r2-3N-2)  +  3(2r2-N-l) 

1  (22) 

+  r(q+2)  (6p^+6r2-3N+r+qr)  +  ij 

By  replacing  r^^  by  r^  and  p^  by  p^  in  equation  21,  we  obtain 
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-  [^2(p2+1)  (2p2+6r2-3N-2)  +  3(2r2-N-l) 
r(q+2)  (3p2+6r2-3N+r+qr)  +  ij 


The  preceding  results  are  also  applicable  to  D^-,  D-_,  D„,,  and 
The  resulting  expressions  are 


‘^11’  °12’  °16’  °Z2’  *^26’  °66^  "  3  |  V  J 

*  lSa\.  (V,’  <VJ'b  (24) 

•  f(Q„)  .(V  ,0,  (5jj)  .0.  (^  )  ICfl 

L  90  90  ^  90  90j  1 


Equation  24  represents  the  closed-form  sums  for  the  flexural  stiffnesses 
of  a  laminate  con^josed  of  q  repeated  sets  of  [0  /±$  /90  ]  laid  up 

^1  ^2 


s>’iiTnsotricaHy  relative  to  the  midsurface  of  the  laminate.  These  stiffnesses 
are  required  to  conduct  analyses  of  orthotropic  laminated  plates.  The 
equivalent  extensional  constants  are  obtained  by  combining  equations  13  and 
14  in  the  appropriate  manner. 

For  cases  when  a  more  general  laminate  layup  is  under  consideration, 
equation  9,  which  is  valid  for  any  laminate  ply  orientations,  should  be 
utilized. 


CALCULATIONS  OF  EXTENSIONAL  AND  FLEXURAL  CONSTAhlTS 


A  computer  program,  AC- 2,  has  been  developed  to  perform  the  calculations 
required  to  determine  the  extensional  elastic  constants  and  flexural  stiff¬ 
nesses  for  a  balanced  laminate  laid  up  symmetrically  relative  to  its  mid¬ 
surface.  This  program  incorporates  the  closed-form  solution  for  D. .'s 
described  by  equation  24  and,  as  such,  is  ver>'  rapid  for  thick  as 
well  as  thin  laminates.  The  elastic  constants,  as  defined  by  equation  13, 
for  several  laminate  families  of  Narmco  5505  as  obtained  from  AC-2  are  plotted 
in  figures  6  through  17  for  room  tenperature  (RT)  and  350°  The  elastic 
constants  of  eadi  lamina  in  its  natural  coordinate  system  (L,  T,  Z)  were 
obtained  from  reference  2. 

The  elastic  constants  that  can  be  found  from  these  curves  are  Ex,  Ey, 

Gxy,  ^yx*  To  assist  the  designer  further  in  choosing  the  desired 

laminate  orientation,  lines  were  drawn  which  represent  both  "I  at  QO"  and 
at  900,"  even  though  only  one  need  be  given  to  derive  the  other. 

Also,  it  should  be  noted  that  the  Gxy  curv^es  are  independent  of  the  "%  at 
Oo"  versus  the  "%  at  90O"  plies.  That  is,  as  long  as  the  "I  at  ±0O'-  values 
are  known,  then  the  remaining  relative  percentage  of  the  0°  versus  the  90O 
plies  does  not  affect  Gxy.  Furthermore,  considering  the  number  of  plies  at 
OO,  ±00^  and  90O  to  be  nl,  n2,  and  n3,  respectively,  then  it  is  seen  that 

1.  No  curve  is  given  for  Ey  for  [0jji/i45n2/90n3]  laminate  orientations. 
This  is  because  Ey  for  these  orientations  can  be  found  by  using  the 
c.ur\^es  for  Ex  after  interchanging  the  at  0®"  and  "%  at  90O" 
labels. 

2.  The  curves  for  [0nl/±30n2/90n3]  laminate  orientations  are  not  given, 
but  the  curves  for  the  [0nl/±60n2/90n3]  laminate  orientations  can  be 
used  for  f0nl/±30n2/90n3]  properties  if  the  following  are  considered: 

a.  Ex  for  [0ni/±30n2/90n3]  is  found  by  interchanging  the  "%  at  QO" 
and  "%  at  90O"  labels  and  using  the  Ey  curves  for 
[0ni/±60n2/90n3]- 
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b.  Ey  for  [0ni/±30n2/90n3]  is  fovind  by  interchanging  the  "I  at  OO" 
and  at  90°"  labels  and  using  the  Ex  curves  for 
[0nl/±60n2/90n3]. 


I'xy  for  [0ni/±30n2/90n3]  is  found  by  interdianging  the  "I  at  0°" 

and  "I  at  90°"  l^els  and  multiplying  the  J'xv  value  from  the 

[  0ni/±60j^2/90n3 1  curves  by  the  ratio  E  /E  for 

X  y 


The  program  also  confutes  the  flexural  stiffnesses  of  the  laminate 
assuming  it  to  be  a  homogeneous  orthotropic  medium  and  thus  utilizing  the 
extensional  constants,  e.g.,  =  Ex  tVl2(l-*^  ^'yx)*  sud  conpares  these 

stiffnesses  to  those  obtained  by  the  rigorous  method,  i.e.,  equation  14. 

ITiis  output  permits  both  the  determination  of  how  thick  a  laminate  of  a  cer¬ 
tain  orientation  must  be  before  the  stiffnesses  Di6  and  D26  approach  zero, 
and  determination  of  the  error  associated  with  the  use  of  flexural  stiff¬ 
nesses  Dll,  Di2,  D22»  and  055  based  on  the  extensional  constants  rather  than 
equation  14. 

The  latter  data  permit  an  assessment  of  when  the  classical  solutions  for 
the  response  behavior  of  homogeneous  orthotropic  plate  and  shell  elements  are 
valid  (for  the  orientation  in  question),  i.e.,  cases  where  the  ratios  D16/D11 
and  D26/D11  approach  zero.  A  further  discussion  of  the  significance  of  the 
stiffiiess  terms  Dff)  and  D26  and  the  ratios  Dig/Dn  and  D26/D11  in  obtaining 
static  and  dynamic  response  solutions  to  solid  laminated  plates  is  contained 
in  section  III. 
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%  at  ±45° 

Figure  7.  Laminate  vs  Percent  Laminae  at  o°,  +45°,  and  90°  for  Narmco  5505  at  350 


R,T.  Single  Ply  350“?  Single  Ply 

Unifiirectional  Properties  ^  Unidirectional  Properties 


%  at  ±li5° 

Figure  8.  Laminate  G^y  vs  Percent  Laminae  at  0°,  ±45®,  and  90®  for  Narmco  5505  at  RT  and  350 


%  at  ±60° 

ure  11.  Laminate  Ej,  vs  Percent  Laminae  at  0®,  +60°,  and  90°  for  Narmco  5505  at  RT 


3^0“?  Single  Ply 


Figure  12.  Laminate  vs  Percent  Laminae  at  0°,  tn0°,  and  90°  for  Narmco  5505  at  350 


Laminate  Ey.  vs  Percent  Laminae  at  0°,  +60°,  and  90°  for  N'armco 
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Laminate  vs  Percent  Laminae  at  0°,  ±60°,  and  90°  vor  Narmco  5505  at  350 


Laminate  Gv>r  vs  Percent  Laminae  at  0° ,  +60° ,  and  90°  for  Narmco  5505  at  RT  and  350 


SECTION  III 


PUT  PUTES 


EQUATIONS  OF  MOTION 

The  equations  of  motion  for  the  plate  element  shown  in  figure  18  are 
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where  the  plate  rotary  inertia  is  considered  to  be  negligible,  P  is  the  mass 
density  of  the  plate,  t  is  the  plate  thickness,  the  force  and  moment  compo¬ 
nents  are  defined  by  equation  6,  and  the  condition  which  each  equation 
enforces  for  the  element  is  indicated  within  the  parentheses.  Also,  all  plies 
are  considered  to  have  the  same  density,  which  is  generally  the  case.  TUie 
shear  forces  and  Qy  can  be  eliminated  from  these  expressions  by  combining 
equations  25c,  25d,  and  25e  to  yield  the  relation 


a^i 


2  y  ^..2  xy  dxdy 


dxT  ^  dy 


(26) 


+  2 


dxdv 


dy^ 


at2 


29 


Nx 


N 


5 


+  •— i  dy 


M  y  9y 

I  -/iX'' 

y'  /  fn"^-'<^''^'-l  *' 

-i./^‘  ^  yx  8Y  ' 

^X^TT^dx  8Nyx^  ' 

M  + - '-m.Y  xy  8x  m  dy.,->^ 

X  8x  X  ay 

8N^ 

^x  ^  8x 


Figure  18.  Plate  Element  Subjected  to  Force  and  Moment  Conponents 


I£,  as  before,  the  disc’ission  is  limited  to  a  laminate  with  a  ply  layup 
idiich  is  symmetric  relative  to  its  midsurface,  then  combining  equations  7b  and 
25e  and  neglecting  in-plane  inertia  terms  yields 
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This  equation  is  based  on  the  assunption  that  the  normal  deflections  of  the 
plate  are  small  conpared  to  its  thickness,  and  the  laminate  ply  orientation 
is  symmetric  relative  to  its  midsurface.  The  latter  assumption  permits  the 
uncoupling  of  the  bending  and  membrane  force  and  displacement  conponents. 

An  examination  of  equation  27c  indicates  that  except  for  cases  when  the 
laminate  layup  is  conposed  totally  of  0-  and  90-degree  plies,  the  and  D26 
terms  are  nonzero.  Hoivever,  for  thick  laminates  composed  of  many  repeating 
sets  of  plies,  the  and  D26  stiffness  terms  tend  toivard  zero. 

This  can  be  seen  in  figures  19,  20,  21,  and  22,  which  plot  the  ratio  of 
these  stiffnesses  to  the  and  D22  stiffnesses  for  50 -percent  volume  frac¬ 
tion  boron/epoxy  laminates.  If  we  consider  the  average  ply  thickness  to  be 
0.005  inch,  we  see  that  for  the  laMnate  orientations  considered,  the  ratio 
[Dig/Diil  or  |D26/Di]^|  is  less  than  7  percent  of  Du  for  thicknesses  of  0.20 
inch  (40  plies)  or  greater.  In  addition,  for  th’ck  laminates,  the  flexural 
stiffnesses  can  be  defined  in  terms  of  the  extensional  constants  for  the  basic 
ply  set  \diich  is  repeated  to  form  the  laminate.  This  is  because  each  basic 
set  of  thick  laminates  undergoes  nearly  a  uniform  strain  field  during  bending, 
as  shoivn  in  figure  23,  while  this  is  not  true  for  thin  laminates. 


THIN  LAMINATE  THICK  LAMINATE 


Figure  19.  Thin  and  Thick  Laminate  Deformation 


Number  of  Plies  for  Several  Laminate  Orientations 


Number  o£  Plies  for  Several  Laminate  Orientations 
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Figure  23.  D16/D22  vs  Number  of  Plies  for  Several  Laminate  Orientations 


Thus,  for  thick  laminates  and  zero  temperature  change,  the  nonzero 
flexural  stiffness  can  be  defined  by 


11  12  (1-y  V  ) 

xy  yx"^ 


22  12  (1-1/  p  ) 

xy  yx 


v  L  t 

D  =  D  =  _ 

12  21  12  (1-1/  p  ) 

xy  yx 


G  t‘ 

D  =  - 
66  12 


where  the  elastic  constants  for  each  laminate  are  obtained  from  equation  13. 
Equation  27c  reduces  to  the  following  classical  governing  equation  for  spe¬ 
cially  orthotropic  plates 


ftcV  ^ay^  at2 


+  N,  ^  +  N  ^  +  2  N 


xy  dxdy 


where 


h  -  “ii 


■’s  '  “l2  "  2  “66 


“2  “22 


A  broad  class  of  closed- form  orthotropic  plate  solutions  to  equation  29 
are  available  in  the  literature,  and  are  valid  for  the  analysis  of  thick, 
laminated,  filamentary  con^osite  structures. 
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PRESSURE  PANELS 


The  classical  theoiy  (reference  11)  for  the  deflection  and  moments  of  a 
uniformly  loaded,  sinply  supported,  specially  orthotropic  laminated  plate  can 
be  obtained  from  equation  29  by  setting 


N=N=N  =^=0 
X  y  xy  ^^2 


and  the  plate's  normal  deflection  W  in  the  form 


E  E 

m=l,3,5  n=l,3,5 


.  .  mffx  .  nrr\ 

A  sin - sin— 

mn  a  b 


and  the  uniform  normal  pressure  by 


“EE; 

nFl,3,5  n=l,3,5  ^ 


16  q  .  mTTx  .  nny 

■= — ^  sin - sin  -r^ 

2  a  b 

mn 


then  substituting  equations  30  and  31  into  equation  29  and  solving  for  Aj^n* 
This  procedure  results  in  the  following  relationship  for  the  plate's  normal 
deflection: 


"'16^  E  E 


^  m=l,3,5  n=l,3,5 


where 


16  q  b _ 


The  stresses  of  any  i^'^  ply  at  a  distance  z  from  the  midsurface  of  the 
laminate,  which  correspond  to  the  deflection  field  described  by  equation  32, 
can  be  obtained  from  equation  5  as 


E  E 

^  1  -  in=l,3,5  n=l,3,5 


;r2|nL^  y(k)  ^ 


'mn  ]  2  yx  ,2 


xy  yx 


.  mjTx  .  nffy 

sin - sin  -r — 

a  b 
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(i)  _ 


-  z  ^  o  f  2  . .  2’ 

1  „Ci)  (i)'  4^  r  4-i  r  nin’‘  ,2  xy  2 

1  -  l>'-  ■'  If''  ■'  in=l,3,5  n=l,3,5  -  lb  a  . 

xy  yx 


.  HiTTx  .  n;ry 

sin - sin  -r-^ 

a  b 


£or_0-  and  £0- degree  plies.  For  plies  o£  other  orientations,  the  contributions 
o£  Q25  and  must  be  considered. 

For  the  case  when  in -plane  conpression  £orc^  and  N  are  also  present, 

A  becomes 
mn 


TT  mn  j  D 


-(a'kiff 


16  q  b' 


+  N  n 

y 


/mnb\  ^  _  4 


UNIAXIAL  AND  BIAXIAL  COMPRESSION 

The  classical  solution  £or  the  buckling  of  specially  orthotropic,  simply 
supported  plates  is  obtained  from  equation  29  by  setting 


N  .,=§4=0 


and  taking  w  in  the  form 

,  .  nOTx  .  nTry  /--z^ 

w  =  A  sin - sin 

mn  a  b 

where  m  and  n  denote  the  number  of  half-wave  lengths  along  the  X  and  Y  axes , 
respectively. 

Noting  that  the  forces  N  and  N  are  conpression,  we  obtain  from 
equation  29  ^  Y 

k  77  V. 

CN  )  = 

^•cr  b^ 


. . 


where 


^  N 

X 

and  m  and  n  are  the  half-wave  length  integers  which  miniirdze  k^*  For  the 
case  when  only  is  present,  i.e.,  <f>  =  0,  and  the  plate  is  very  long 
(a/b  »  1) ,  we  have,  after  noting  that  kx  is  minimized  in  this  case  for  n  =  1, 


or  (k  )  occurs  when 

X  . 
rain 


when  m  is  a  whole  number  and  we  see  for  a  long  plate  that  k^  becomes 


k  =  2 

X 


1  + 


'”12*2  V 


V 


”n  ”22 


(38) 


Figures  24  and  25  are  plots  of  Nxcr  for  [0/±45/0]2s  laininates  at  roran 
tenperature  and  350°  F  as  a  function  of  panel  aspect  ratio,  panel  width,  and 
ratio  of  Ny  to  Nx-  These  curves  were  obtained  from  equation  37. 


b  (Inches) 

Figure  24.  Buckling  Load  Nx,cr  vs  Width  b  for  a  Siinply  Supported  Plate  With  My  (RT) 


natural  FREQIJENCIF.!=; 


■node  f^-Kncies  associated  with  various  m  and  n 

mode  shapes  for  a  singly  supported  plate,  w  is  taken  in  the  form 


w(x,  y,  t)  =  A  sin  ^^^sin  ^ 

mn  a  b 

where  6>  is  the  angular  velocity  in  radians  per  second.  Substitution  of 


(39) 


(40) 
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SECTION  IV 


ANALYSIS  OF  ORIHOTROPIC  SANDWICH  PLATES  AND  SHELL  STRUCTURES 

CONSTITUTIVE  RELATIONS 


In  solid  plates  and  shells,  shear  strain  is  generally  neglected.  But  in 
sandwich. plates  and  shells,  shear  strain  in  core  needs  to  be  considered.  The 
sandwich  element  is  conpised  of  two  face  plates  and  a  core  as  shown  in  fig¬ 
ure  26.  The  distance  between  the  centroids  of  the  two  faces  is 

h  =  c  +  t^  (41) 


Consider  a  cylindrical  sancfttfich  shell  element  as  shoi>m  in  figure  27.  Let 
X,  Y,  Z  be  a  set  of  orthogonal  axes  itfith  X  along  the  genatrix  of  the  cylinder, 
Y  tangent  to  the  cylindrical  surface,  and  Z  normal  to  the  cylindrical  surface, 
and  u,  V,  w  be  displacements  along  these  three  axes,  respectively.  Let  N^^,  N^, 
Njjy,  denote  the  sectional  forces;  My,  Mxy  the  sectional  moments;  and  4c, 

Qy,  the  transverse  shear  forces  per  unit  length  of  the  element.  This  element 
is  subject  to  a  distributed  lateral  load  q  per  unit  area. 

The  equations  of  motion  of  this  element  are  witten  as 


dM  dM 

_y  +_J2L.  n 

5x  6y  X 


=  I 


(■12) 


aM  aM 

ay  ax  y 


=  I 


at' 


(43) 


N  2N 

•’'5x2 


xy  axoy  R  ax  ay 


=  M 


a^w 

at^ 


(44) 


where  I  denotes  the  rotary  moment  of  inertia,  M  the  mass  per  unit  area  of  the 
element,  R  the  radius,  t  the  time,  and/Jy  rotations  as  shown  in  figure  28. 


Let  and  denote  the  mass  densities  of  the  face  and  the  core,  respectively. 


I 


(45) 


M  =  2t^P^.cp; 


(46) 
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NATURAL  FREQUENCIES 


If  we  seek  the  fundamental  frequencies  associated  with  various  m  and  n 
mode  shapes  for  a  simply  supported  plate,  w  is  taken  in  the  form 


r  ,  •  m^x  .  nTry 

w(x,  y,  t)  =  e  ""A  sin - sin 

mn  a  b 


where  is  the  angular  velocity  in  radians  per  second.  Substitution  of 
equation  39  into  equation  29  and  setting  the  in-plane  forces  and  normal 
pressure  q  to  zero  yields 


I 

1 

I 
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SECTION  IV 


ANALYSIS  OF  ORTHOTROPIC  SANDWICH  PLATES  AND  SHELL  STRUCTURES 

CONSTITUTIVE  RELATIONS 


In  solid  plates  and  shells,  shear  strain  is  generally  neglected.  But  in 
sandwich,  plates  and  shells,  shear  strain  in  core  needs  to  be  considered.  The 
sandwich  element  is  composed  of  two  face  plates  and  a  core  as  shown  in  fig¬ 
ure  26.  The  distance  between  tho  centroids  of  the  two  faces  is 

h  = 


Consider  a  cylindrical  sandt/ich  shell  element  as  sho\^nl  in  figure  27.  Let 
X,  Y,  Z  be  a  set  of  orthogonal  axes  with  X  along  the  genatrix  of  the  cylinder, 
Y  tangent  to  the  cylindrical  surface,  und  Z  normal  to  the  cylindrical  surface, 
and  u,  V,  w  be  displacements  along  these  three  axes,  respectively.  Let  Nj^,  N^., 
Nj{y,  denote  the  sectional  forces;  My,  M^y  the  sectional  moments;  and  4c, 

Qy,  the  transverse  shear  forces  per  unit  length  of  the  element.  This  element 
is  subject  to  a  distributed  lateral  load  q  per  unit  area. 

The  equations  of  motion  of  this  element  are  witten  as 


aM  aM 

ax  ay  X 
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(42) 
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xy  ax6y  R 
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^2 

a  w 
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(44) 


where  I  denotes  the  rotary  moment  of  inertia,  M  the  mass  per  unit  area  of  the 
element,  R  the  radius,  t  the  time,  andySy  rotations  as  shown  in  figure  28. 


Let  and  p'f-  denote  the  mass  densities  of  the  face  and  the  core,  respectively. 


I 


(45) 


M  =  2tjP^  +  cp^ 


(46) 
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Referring  to  figure  28,  we  see  that  the  u  and  v  components  can  be  taken  in 
the  form 

u(z)  =  Bz\  vCz)  =  Bz  (47; 

X-  y 

Assume  the  bending  moment  to  be  taken  solely  by  the  faces  and  shear  solely  by 
the  core,  then 


^v* 

«  g' 

/8v.3w\ 

=  n' 
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^cy 

\^"ay; 

^cy 

(  y  ay) 

^xz 

=  g' 

cx 

/8u  8w\ 
iSz  Sxj 

=  g' 

cx 

The  face  sheets  are  considered  to  be  specially  orthotropic  and  thus  their 
stress-strain  relation  is  given  by 


a  =  T - (e  +v  e  ) 

XX  1-v  V  ^  X  yx  y 
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Similarly, 


5  ®i2L) 
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The  section  raranents  are  obtained  by  integration. 
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where 


"'ll  " 
Similarly, 
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Substituting  equations  53,  54,  and  55  into  equations  42  and  43  we  obtain 
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V0x‘ 


yxa^ay , 
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Now  let  us  consider  the  rase  of  a  simply  supported  rectangular  plate. 


Let 
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From  equations  49,  58,  59,  and  60,  we  find 
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Substituting  equations  59,  62,  and  63  into  equation  56  and  simplifying  yields 
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Equation  64  can  be  written  in  the  following  compacted  form: 

a,T  A  +  a,-  B  +  a,  -  C _  =  0 

11  mn  12  mn  13  mn 


(64) 


(65) 


(66) 


where  a,,,  a,-,  and  a,_  are  the  coefficient  O’  A__,  B  ,  and  C  .respectively. 
11  12  13  '  mn  mn  mn 


Similarly,  substituting  equations  59,  62,  and  63  into  e^^uation  57,  we 
obtain 
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where 
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A  plate  can  be  considered  as  a  circular  shell  with  infinite  radius .  R  in 
equation  44  becones  infinity.  Let  1^  =  0.  Then  eliminating  Qx  and  m 
equation  44  by  equations  56  and  57,  we  obtain 


Substituting  equations  58,  62,  and  63  into  the  preceding  equation  yields 
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Fr<Mn  equations  66,  67,  and  70,  we  obtain 
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Let  us  define  the  following  nondimensional  ratios: 
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Then  the  coefficients  of  equations  66,  67,  and  70  may  be  written  as 
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COMPRESSION  STABILITy 

When  the  loads  and  Ny  reach  critical  levels,  the  denominator  of  equa¬ 
tion  74  vanishes.  This  yields 


a„  +  mSl  +  h^N  =  (77) 

31  X  y 

■(^33 *-^12^21  "  ^11^22^  ^32^®11®23'®13^21^) 

*  -  -  =  -a^j 

^12^23  ■  ®22®13 

Let  Njj  and  Ny  be  positive  for  caiqiression;  the  signs  of  Nx,  My,  are  changed 
in  the  preceding  equation,  giving 


_2  2 
in  N  +  h  N  = 
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From  the  last  two  eqitations  of  equation  75  and  the  preceding,  we  have 
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Hie  critical  buckling  load  is  obtained  (equations  80,  81,  and  82)  by  finding 
the  integers  m  and  n  idiich  give  mininuB  values  of  kjj  or  ky.  These  values  for 
the  cases  of  uniaxial  and  biaxial  conpressiwi  are  shnm  in  figures  29 
through  36. 
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Figure  29.  Buckling  Coefficient  Kj^ 


^  ^  ?T^tf/e)(mx/c) 

2  G^Wc)  * 


V  for  Simply  Supported  Flat  Sandwich 
Panels  (RT) 


(Nx)cr 


y  _  7T^’A(Vc)(3.nj./c) 
2  G^(D/c)  2 


Figure  31,  Buckling  Coefficients  for  Biaxially  Loaded  Simply  Supported  Flat 

Sandwich  Panels  of  a/b  =  1  (RTl 
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2  G^CVc)^ 


0  .  .5  1 


K 

y 

Figure  33.  Buckling  Coefficients  for  Biaxially  Loaded  Siii?)ly  Supported  Flat 

Sandwich  Panels  of  a/b  *  3  (RT) 
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Figure  34,  Buckling  Coefficients  for  Biaxially  Loaded  Simply  Supported  Flat 

Sandwich  Panels  "of  a/b  »  1  (350“  F) 


Figure  35.  Buckling  Coefficiaits  for  Biaxially  goaded  Singly  Supported  Flat 

Sandwich  Panels  of  a/b  ■  2  (350“  F) 


2b^ 

msma 


Figure  36.  Buckling  Coefficients  for  Biaxially  Loaded  Siji5)ly  Supported  Flat 

Sandwich  Panels  of  a/b  =  3  (350°  F) 


In  a  similar  manner,  the  natural  frequencies  o)  for  any  m,  n  mode  shape 
can  be  obtained  by  solving  equation  78  for  given  and  ky.  As  pointed  out  by 
previous  investigators  (reference  ^ ) ,  there  are  three  resonant  frequencies 
for  each  mode  instead  of  a  single  resonant  frequency  for  solid  plate.  TWo  of 
these  resonances  are  associated  with  the  two  rotations ,  and^y,  and  one  is 
associated  with  lateral  deflection  6>.  For  the  cases  where  the  core  is 
isotropic  (G^x  =  the  three  resonant  frequencies  reduce  to  two  for 

each  mode  of  vibration. 


For  filament-reinforced  face  sheets,  the  following  expressions  are 
invariant  with  respect  to  the  rotation  of  the  reference  axes  in  the  plane  of 
the  plate  (reference  4). 


l-y  p 
xy  yx 


1-1/  V 


where  L  and  'f  denote  the  longitudinal  and  transverse  directions  of  the 
filaments . 


Since  the  (i,j  =  1,  2,  3)  depend  on  the  five  face-sheet  elastic  constants 

through  Dll,  ^22\^66>  the  critical  values  of  and  Ny  depend 
on  these  five  variables.  The  buckling  constants  k^  and  ky  are  the  ratios  of 
Nx  and  divided  by 


and  hence  are  defined  by  four  independent  variables.  (Only  four  of  the  five 
elastic  constants  E^,  *'xy»  *yx»  Scy  independent  because  of  the 

auxiliary  relationship  Ex^yx  using  the  invariant  relations 

(equations  83  and  84),  k^  and  ky  may  be  expressed  in  terms  of  Ii,  I2,  and  two 
variables .  Let 


We  can  express  y  and  P  in  terms  of  a  and  ^  as 


(85) 


(86) 


6Q 


T 


1 


¥ 


IS 

% 


■s' 


- 


2  {'^ 


-a 


1.1 


f87) 


For  a  given  composite  face  sheet  on  a  given  core  with  given  c,  and 
the  variation  of  critical  kj^  and  ky  may  be  plotted  against  the  properties  of 
the  core  for  different  values  of  a  and  tU.  This  reduces  greatly  the  computa¬ 
tion  and  presentation  of  the  critical  values  of  k^^  and  ky  for  different 
orientations  of  the  face  plates.  ^ 


A  coiH)uter  program,  AC-5,  entitled  "Honeycomb  Panel  Stability  Under 
In-plane  Biaxial  Loading,"  has  been  developed  to  perform  the  calculations 
indicated  by  equations  80  and  81.  This  program  searches  and  finds  the  half- 
wave-  length  numbers  m  and  n  which  minimize  the  buckling  load  for  a  given 
panel  geometry  and  face  sheet  material.  The  .^C-5  program  was  utilized  to 
generate  the  stability  curves  for  the  [02/i45]^  laminate  orientations  used 
for  the  structural  elements  tested  in  this  program  and  reported  in  volume  IT. 

For  the  Narmco  5505  system,  the  values  of  the  material  invariants  used 

were 
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00 

00 

Msi 

■2 

= 

0. 

128 

Msi 

= 

31 

.54 

Msi 

= 

0. 

066 

Msi 

Figures 

29 

and 

30 

give 

for  room  temperature 


for  350°  F 


and  350°  F,  respectively,  for  the  [02/±45](^  orientation.  Figures  3]  through 
36  give  'the  buckling  loads  for  the  same  orientation  at  room  temperature  and 
350°  F  for  various  combinations  of  longitudinal  and  transverse  in-plane  loads 
and  aspect  ratios. 


UNIFORMLY  DISTRIBUTED  PRESSURE  .AND  EIXi'E  CO^IPRESSIOX 

We  now  consider  that  the  previously  discussed  panel  is  subjected  to  a 
uniform  lateral  load  q  as  shoun  in  figure  37.  The  pressure  q  can  be  expanded 
into  the  following  Fourier  sine  series: 


q 


EE 

m  n 


.  m;rx  .  n;ry 
Qmr,  Sin  —  sin 
^mn  a  b 


(88j 


% 
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where 


=  16q/7r  mn 
%n  "  ^ 


for  odd  values  of  m  and  n 
for  m  or  n  even 


y'  "x 

X 

Figure  37.  Sandwich  Panel  Subjected  to  Normal  Load  and  Edge  Ccnqjression 


In  this  case,  the  loading  is  static  and  the  radial  velocity  is  zero,  i.e., 
w=  0;  thus,  the  panel  deflection  is  expressed  as 


CC  X 

E  E 

iiv=l,3  n=l,5,5 


,  .  mjrx  .  uTTv 

A  sin - sin  -r-' 

mn  a  b 


and  we  obtain  from  equations  74  and  88  that  takes  the  form,  with  and 
Ny  taken  to  be  positive  for  compression: 

16  q  b^  A 

A  =  - (S 


where 


(®12^23  ■  ^22^13) 


+R  k  n 
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In  a  similar  manner*  we  find  that 
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dy 


TT^  in=l,3  n=l,3 


CO 
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.  -I  1  mn 


,  mb 


V  ^  ^  rV  (b/_al  _  ^  (^'j 
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1 
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I  b  ,,  2  b 

C^-2y)IX^  m-V  -  m  — 


Sin  cos  ^ 

(100) 

,ere  the  Eoltaing  auxiliary  relationships  were  used  to  reduce  the  preceding 
::juations  to  their  final  form. 


E.  9 

y  =1/  ^  ^  V  a 

"xy  yx  E  yx 


(101) 


^yx“=^-^’' 

The  stresses  in  the  panel  are  then  detennined  from  the  relationship 
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’■xt-V'^-’-yt'V'^ 


^x""t^(c+t^)  2t^ 
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(102) 
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nd,  as  mentioned  before,  the  in -plane  forces  Ny  and  N,  are  taken  to  be 

jositive  if  they  are  compressive. 

trr-  n,  Af  7  has  been  developed  to  perform  the  calculations 
A  computer  S  internal  moments  and  shears  of  simply 

to  obtain  the  noimal  deflecuons  md  intero  ^  subjected 

supported  honeycomb  TabUs  I  and  II  contain  data  obtained 

KrarruS^Iertfiested  in  this  ptogtam  and  teported  in 
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volume  II  of  this  report.  Figure  38  is  a  plot  of  the  maximum  normal 
deflection  for  panels  with  face  sheets  composed  of  [02/±45]q  laminates  and 
no  in-plane  forces.  These  data  were  obtained  from  AC- 7. 


TABLE  I.  INTERNAL  MOMENT  AND  SHEAR  COEFFICIENTS  FOR  SIMPLY 
SUPPORTED  PANELS  WITO  NO  IN-PLANE  LOADS* 


a/b 

Temperature 

M  /qb^ 

x/a  =  0.5 
y/b  =  0.5 

M  /d? 

y 

x/a  =0.5 
y/b  =0.5 

Q^/qb 

x/a  =  0.0 
y/b  =  0.5 

Qy/qb 

x/a  =0.5 
y/I)  =  0.0 

1.0 

RT 

0.073 

0.027 

0.396 

0.262 

3500  F 

0.072 

0.026 

0.396 

0.254  ! 

2.0 

RT 

0.094 

0.077 

0.485 

0.408 

3500  F 

0.104 

0.071 

0.395 

3.0 

RT 

0.092 

0.104 

IHl 

0.461 

3500  F 

0.107 

0.099 

U.53.S 

0.451 

*  Laminate  orientation  -  [02/±45]^ 


Material  -  Narmco  5505 


G' 

cx 

G' 

cy 
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Figure  38.  Parameter  w  ^/qb'^  vs  V  for  Simply  Supported  Honeycomb  Sandwich  Panels 


68 


I  m  ‘  « 'll  ' '  •"  ’5’  ■  I 


TABLE  II.  MAXIMUM  NORMAL  DEFLECTION  AND  INTERNAL  MOMENTS  AND  SHEARS 
FOR  SIMPLY  SUPPORTED  PANELS*  WITH  N  /(N  )  =0.333 


X  X  cr 


a/b 

Temper¬ 

ature 

x/a  =0.5 
y/b  =0.5 

M  /qb^ 

x/a  =0.5 
y/b  =0.5 

M  /qb^ 

y 

x/a  =0.5 
y/b  =  0.5 

Q^/qb 

x/a  =0.0 
y/b  =0.5 

x/a  =0.5 
y/b  =  0.0 

1.0 

RT 

0.00581 

0.107 

0.044 

0.562 

0.356 

3500  F 

0.00547 

0.111 

0.040 

0.573 

0.342 

2.0 

RT 

0,0290 

0.121 

0.149 

0.850 

0.594 

3500  F 

0.0308 

0.148 

0.145 

0.864 

0.597 

*Laminate  orientation  -  [02/±45]^ 
Material  -  Narmco  5505 

C  “  DQv  “ 


■  DQy  " 

^  “  y/^ll  022 

SHEAR  STABILITY  OF  FLAT  SANDtVICH  PANELS 

A  computerized  method  has  been  developed  to  predict  the  shear  stability 
allowable  o£  flat  sandwich  panels  with  orthotropic  face  sheets  and  core.  The 
face  sheets  are  assumed  to  be  similar,  and  the  shear  loading  is  in  the  plane 
of  the  panel.  The  analytical  developnent  is  given  in  reference  5,  and  only 
the  final  governing  equations  are  presented,  along  with  the  listing  of  the 
computer  program  and  a  sample  problem  solution. 

Consider  the  sandwich  panel  to  lie  in  the  X-Y  plane  with  the  Z-axis 
perpendicular  to  the  face  sheets,  as  shown  in  figure  39.  The  critical  shear 
buckling  load  is  given  by 


where 


CNxyVr  = 


l*'^xy''yx 


(C+tf)2 


(103) 


and  the  other  terns  are  defined  in  the  figure. 
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T 


T 


*Jcx>  ^cy  ~  core  shear  moduli  in  the  x  and  y  direction,  respectively. 
Using  the  contracted  notation  that 


i  =  n  +  q(m  -  1) 
j  =  s  +  q(r  -  1) 

where  q  is  the  limit  of  m,  n,  r,  and  s,  the  resulting  governing  equation 
becomes 


where 


KA.C.  +  y  H..C. 
1  1  Z-rf  ij  j 

j=l 


■"I 


=  0 


Dividing  through  by  and  rearranging  terms  results  in 


q  H. .C. 
K  c.  =  -  V'  -^J  J 


-i  ■  ■  £  A. 
j=l  ' 

whicli  in  matrix  notation  becomes 

b]  [c|  -  k|c||  =  0 


(105) 


(106) 
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It  is  noticed  that  the  last  equation  is  in  the  form  of  tlic  standard 

eigenvalue  problem.  For  this  typo  of  problem,  a  modified  matrix  iteration 

process  (reference  7)  will  yield  the  highest  eigenvalue,  which  mc.'ins  tliat  we 

have  the  latest  value  of  K  since 

s 


To  facilitate  the  solution  for  K,  it  is  pointed  out  in  reference  5  tiiat 
the  matrix  [B]  can  be  separated  into  oven  and  odd  terms  of  Cjiin>  i-c.,  [1^1  even 
for  m  +  n  even  and  [B]  odd  for  m  +  n  odd.  Then,  [B]  even  and  [I^]  odd  can  cacli 
be  solved  independently  for  Kg  even  and  Kg  odd,  and  the  minimum  of  these  two 
values  would  be  the  same  value.  Kg,  that  would  result  if  the  original  iB) 
matrix  were  utilized  to  find  Ks.  This  procedure  greatly  reduces  computing 
time. 
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SECTION  V 


STIFFENED  SKIN  CONSTRUCTION 


CONSTITUrn^  REUTIONSUIPS 

The  analysis  in  this  section  for  an  K-ply  flat  sheet  stiffened  with  Z  or 
hat  stiffeners  parallels  the  work  performed  by  Lacknan  and  Ault  (reference  6) , 
but  is  extended  to  account  for  a  heterogeneous  sheet  rather  than  tlie  homogeneous 
sheet  treated  by  these  authors.  We  noitf  seek  to  analyze  the  t>'pical  stiffened 
sheet  shown  in  figure  40  ^dien  the  coordinate  origin  is  located  in  the  midplane 
of  the  sheet. 


Figure  40.  Typical  Stiffened  Sheet 


From  equation  5  we  have  the  following  relationship  between  tlie  stress  and 
displacement  con^wnents  of  the  i^ply  of  the  stiffened  sheet  shown  in  figure 
40: 


And  if  we  further  consider  each  sheet  to  be  cwnposed  of  ply  orientations  which 
are  s>Tiinetric  relative  to  its  midsurface,  then  we  have  from  equation  7  the 
folloiving  relationship  betiveen  the  sheet  in-plane  forces: 


\l  \2  \6  'I 


^2  ^22  ^26 


^16  ^26  he 


d\/dy  ^ 
dy  ^  ^ 


where 


13  “  EM  (N.*!'”!) 


and  is  defined  in  equation  2. 

The  axial  force  per  unit  width,  Njj,  along  the  X  axis  of  the  stiffened 
sheet  is  obtained  from  the  relationship 


b  /2  r 
s  -t/2 


f  «  /  a  dz  +  I  <7  d:  ►  dy  =  +  b  N  ^  (108) 

J  J  X  J  X  ’  XXX 


\n  -(h^-t/2)  -t/2 


stiffener 


when  the  and  denote  the  force  per  unit  width  carried  by  face  and 
stiffener,  respectively.  If  we  now  consider  is  of  the  foim 


_  2 

_  If  _ ®  i  If  ^ 

-  ‘‘i  ex  *  ''2  . 2 

dx 


(109) 


then,  by  combining  equations  109  and  7  with  equation  108,  we  obtain,  after 
integrating, 
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s  ox 


(110) 


Noting  that  in  the  Y  direction  only  the  sheet  carries  the  N»,  and  in-plane 
loads,  we  obtain  from  equation  7 


N  =  N 


f  _ 


^2 


ax  ^22  ay  ^ 


/au  av^ 


26  \ay 


ax 


xy 


Nf 

xy 


au 
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26  ^y 


+  A, 


66  \Qy 


ax 


(111) 


j 


'T 


£ 


T 


I£  v\re  further  consider  the  bending  moment  component  in  the  X  direction  carried 
by  the  stiffener  to  be  of  the  form 


au  .2 

X  nax  4-2 

ax 


then  we  find 


b  M 
s  X 


b  /2 
s 


=  / 


-b  /2 
s 


t/2 


t/2 


J'  (T^zdz  +  J'  zc 

[-{"s  "  '■'2)  -t/. 


dy 


stiffener 


=  +  b 

X  S  X 


Thus,  the  bending  moment  components  for  the  stiffened  sheet  become 
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Substitution  of  equation  118  into  114  yields 


M 

X 


? 

g"w 


16  dxdy 


(119) 


If  the  plate  is  now  subjected  to  an  in-plane  compression  force  and  if 
we  neglect  the  changes  in  the  in-plane  forces  due  to  the  plate's  buckling 
displacements,  the  governing  equilibrium  equation,  for  the  plate  shoun  in  fig¬ 
ure  18  is  obtained  from  equation  26  as 


dh\  a^M  2 

_ X  +  2 _ ^  ^  ^ 

^2  axay  ^2  -x  ^^2 


(120) 
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Substitution  of  equations  115,  116,  and  119  into  equation  120  yields 


D  +  4  D  +  2  /D  +  2  D  \ 

lax^  16^5^3^  [  12  66/  ^^2^2 


dx^ 


dx^'ay 


+  4D  -  +  D 

26  3  22  ^  4 

Sx^  dy 


=  -N^ 


X  ^  2 

ox 


where 


c;  =  Cj/b^ 


<^2  = 


(121) 


By  comparing  equations  121  and  27c,  it  can  be  seen  that  the  solution  to  ajiy 
anisotropic,  laminated,  stiffened  plate  can  be  obtained  from  knoim  solutions 
of  equation  27c  simply  by  substituting  Dj  for  For  example,  for  the  case 

of  a  singly  supported  wide-column  construction,  we  take  w  in  the  form 

w  =  A  sin  (122) 


and  we  find  by  siobstituting  the  chosen  form  of  w  into  equation  12  that  the 
critical  buckling  Iqad  is 


(123) 


i.e.,  the  classical  solution  has  been  modified  by  replacing  with  D]^. 


ZEE  STIFFENERS 


^  ^  To  utilize  equation  121  for  Z-stiffened  sheets,  the  values  of  c{, 

C2,  Ki,  and  must  be  established.  Thus,  we  now  consider  the  typical 
Z-stiffener  shown  in  figure  41.  Referring  to  equation  108,  we  see  that 


/  / 

*/ 

-K4) 

-*s  ■(’'s  * 

-t/2 


j  »,dzd)r 


) 


z 


(124) 


With  the  assunption  that  the  stiffener  only  resists  loading  in  its  stiffening 
direction,  we  see  from  equation  5  that 


(125) 


where  E^ 


is  the  modulus  of  the  stiffener  in  its  stiffening  direction. 
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Combining  equations  124  and  125  and  performing  the  indicated  integration 
results  in 


\vhere 


au  A 

_  Y  _ o 

^1  ax  ^2  „  2 
ax 
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In  a  similar  manner,  we  find  that 


Jc^zdzdy 

-4+  h 

2  s 


n  ax  ^-2  -  2 
ax 


where 


C,  =  -E'  /f  t 
1  X  s  s 


(127) 


HAT  STIFFENERS 

In  the  case  of  hat  stiffeners,  we  consider  the  t>'pical  stiffener  sho\>(n 
in  figure  42.  With  the  use  of  figure  42,  we  find  that 
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SECTION  VI 


MINIMJM  I'/EIGHT  .ANALYSIS 


GENERAL 

A  conqjuter  program,  AC- 12,  was  developed  to  calculate  for  hat-  and 
Z-stiffened  sheets,  as  well  as  the  buckling  parameter  NjjL2  of  equation  115 
for  wide  columns. 

For  tfie  case  of  axial  compression,  the  program  equates  the  local  buckling 
loads  of  the  stiffened  plate  elements  to  insure  that  these  elements  fail 
simultaneously.  Then,  by  equating  this  local  panel  strength  to  the  panel 
strength  in  general  instability,  one  may  obtain  a  minimum  'ceight  concept  for 
a  given  load  requirement  according  to  the  simultaneous  failure  mode  h>'pothesi.s 
for  minimum  weight  structures. 

The  local  buckling  allowables  are  obtained  from  tiie  folloicing 
considerations ; 

1.  For  elements  with  three  sides  simply  supported  and  one  of  tlie  unloaded 
sides  free,  die  buckling  load  can  be  obtained  frcrni 

N  =  k_  tVf^  (131) 

X  7  s  s 

for  Z- stiffeners  with  ky  obtained  by  solving  the  transcendental 
equation  presented  in  reference  6. 

2.  For  elements  v.-'ith  all  four  sides  sinply  supported,  we  have  from 
equation  37  that 


X 


kf  zl/bl 


N 

X 


(152) 
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where 


and  the  symbols  s  and  £  denote  stiffener  and  face  sheet,  respectively. 


ZEE  STIFFENERS 


For  simultaneous 


buckling  of  the  Z-sti££ened  elements,  we  require 


k 


7 


(155J 


or 


A 


while  for  simultaneous  buckling  of  the  face  and  stiffener  elements,  we  require 


(134) 


Equating  the  applied  load  in  the  stiffener  to  the  allowable  stiffener 
load  yields 


N  t 
X  X  s 

,  (2f  +  h  \ 

E^t^  +  E®  i — V — ^t 

X  f  X  b  s 

s 


k  t  ' 
s  s 

h' 

s 


(135) 


where  is  the  applied  compressive  force  to  the  panel.  With  the  use  of  these 
relationships,  we  find 


/  N  E®  b  ^  f 

h  =b  ■?  -  t. 


0  “ 


l^x  ‘f 


t  =  c<  h  tVb 
s  1  s  f  s 


(136) 


f  =  h 
s  Is 


For  a  given  applied  force  and  b^,  computer  program  AC-12  utilizes  the  rela¬ 
tions  of  equation  136  to  calculate  (^^^^ allowable’  optimum  configura¬ 


tion  IS 


(N  L  ) 

i-X-J  =  N 
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HAT  STIFFENERS 


In  a  similar  manner,  we  find  for  hat  stiffeners  that  f  =  h  /2  and 

s  s 
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N  E^ 

X  X 


k  «  ‘ 
s  1 


2  tf 


ll/2 


f  = 
s 


h  /2 
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(137) 
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Using  these  relationships,  program  AC-12  calculates  the  desired  properties  for 
the  hat-stiffened  sheet  in  a  manner  identical  to  the  Z-stiffened  sheet. 


Figure  42.  Typical  Hat  Stiffener 
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SECTION  VII 


THERMDELASTIC  RELATIONSHIPS 


In  this  analysis,  it  is  assumed  that  the  elastic  and  thermal  response  of 
each  ply  of  an  N-ply  laminate,  is  kno\^m.  Tne  quantities  specifying  the  elastic 
and  thermal  response  of  a  ply  may  be  obtained  from  the  data  presented  in 
volume  I.  The  thermoelastic  stress-strain  relation  for  the  i^^  layer  in  a 
state  of  plane  stress  (stress  components  normal  to  the  L-T  plane  are  zero) 
is  given  by 


’T  T  0  1 
11  12 

hi  hi " 

j 

AT 

^LT, 

|i 

“  hi 

ii 

where 


hi  =  h/(i  -  "lt  "tl)  '  \  hl/h 

l(^  '  At  Ti.)  °  At  h  ‘  At  Al) 


T  =  f  E 
12  TL  L 


"^33  “  ^LT 


^LT  ^  ^LT 


(138) 


Figure  43  shows  the  laminate  in  question. 

Here,  denotes  the  magnitude  of  the  strain  in  the  longitudinal  direc¬ 
tion  due  to  a  unit  strain  in  the  transverse  direction;  denotes  the 

magnitude  of  the  strain  in  the  transverse  direction  due  to  a  unit  strain  in 
the  longitudinal  direction. 

To  obtain  the  stress-strain  relations  for  a  lamina  for  the  case  where 
the  applied  stress  makes  an  angle  9^  to  the  natural  lamina  axis,  equation  138 
must  be  subjected  to  a  coordinate  transformation.  The  case  under  consideration 
is  illustrated  in  figure  44. 


91 


The  stress-strain  relations  for  the  typical  i^^  ply  of  the  laminate  shown 
in  figure  43  is  assumed  to  be  of  tlie  form 


'xyj  i 


B  e  -  a  AT 

1  y  y 


(139J 


€  -a  ZiT  I  . 
xy  xy  J  1 


'HiH 


where  the  matrix  [BJi  of  order  3  x  3  is  yet  to  be  determined.  It  is  \\'ell 
kno\m  that  the  ply  stress  and  strain  components  transform  from  one  coordinate 
system  (X,  Y,  Z)  to  another  coordinate  system  (L,  T,  Z)  according  to  the 
transformation  law  of  a  tensor  of  order  tovo  (reference  1).  Hence, 


‘^3cy|  i  I  ^LTli 


^xyli 


(140) 
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where 


2  2 

m.  n.  2  m.n.  1 
1  1  1  r 


n.  m.  -2  m.n. 
3  1  11 


-m.n. 
1  1 


2  2 

m.n.  m.  -n. 
11  11 


(141) 


and 


n.  =  sin  ft ,  m.  =  cos  ft. 

1  i’  1  1 


Substituting  equation  140  into  138  and  performing  the  appropriate  matrix 
manipulations,  we  find  that  [B]^  of  equation  139  is 


[b],  =  HJo] 


(142) 


By  performing  the  indicated  multiplication  for  equation  142,  it  is  found  that 
the  elements  of  the  [B]^  matrix  are 

B^l  =  T^/  .  2  (T^2  ^  2  T33)  m^n^  .  T^^n^ 


®12  ~  ^^11  *  ^22  '  "^12  +  n  )  =  B^^ 


®13  "  ^^11  ■  ^12  ■  2^33^^  *  ^^12  '  ^22  ^  ^33^  "  ^31  (143) 


B22  =  -  2  (T^2  .  2  T33)  nV  . 


®23  ""  *-^11  '  ^12  '  ^^33^  ^  *  ^^12  ~  ^22  ^  ^33^  ""  ®32 


®33  *-^11  ^22  '  ^^12  ’  ^33^  ^33  ^  ^ 
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We  will  now  consider  a  laminate  consisting  of  n  plies  with  natural  axes 
making  various  angles  with  respect  to  the  reference  (X,  Y)  axes  of  the 
laminate.  We  impose  the  restriction  that  the  plies  are  arranged  symmetrically 
with  respect  to  the  center  plane  of  the  laminate.  Because  of  this  restriction, 
there  is  no  coupling  between  in-plane  loading  and  out-of-plane  deformations. 

The  laminate  plate  shoivn  in  figure  43  is  assumed  to  be  subjected  to 
applied  stresses  (ir- plane)  and  temperature  changes  (uniform)  which  do  not 
result  in  curvature  of  the  plies.  For  a  laminate  of  thickness  t,  the  stress 
resultants  are  defined  by 


t/2 

r 

t/2 

t/2 

1  <Xy  dz, 

N  =  Ic  dz, 

J 

y  J  y 

xy  J 

-t/2 

-t/2 

-t/ 

^xv  (144) 


Integrating  equation  139  across  tiie  laminate  thickness  and  assuming  the  stress 
components  to  be  constant  across  each  lamina  and  the  strain  components  not  to 
vary  with  z,  we  obtain 


{■'<}  =  WH  {cl 


(145) 


where 


i  1 

=  Z)  ''xj  hi. 


i=l 


a  h. 

XV  .  1 


®jk  "  S  (^jk)i  ^i 


(146) 


;  +  a  Bt  -  +  «  B. , 

X  11  y  12  xy  13 


)i  hi 


n  . 

=  y^la  B  +ot  Bo,  Vli- 

\  X  21  y  22  xy  23  yi  1 


i=l 


n 


(146) 


xy 


\ 

2-/  fsc  hi  *  “y  ®32  “xy  ®33  ): 
i=l  '  ' 


and  denotes  the  thickness  of  the  ith  ply. 


I*  c  c 

The  coefficiuits  of  thermal  expansion  ,  o(^  ,  for  the  n-ply 

laminate  are  obtained  from  equation  145  by  setting  the  applied  in-plane 
resultants  equal  to  zero.  Hence, 


(147) 


Then  the  general  stress-strain  relation  for  an  n- layered  laminate  subjected 
to  extensional  deformation  takes  on  the  form 


jN|  =  [Bj€|-  =  [®] 


(148) 


Therefore,  for  any  cranbination  of  applied  in-plane  loadings  and  uniform 
terperature  changes,  the  corresponding  strain  field  for  the  laminate  is 
obtained  from  equation  148. 


for  i 


1,  2,  --N  (149) 
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To  calculate  the  stresses  in  the  ply  due  to  the  strain  field  of  equation 
149  for  all  plies,  equation  138  is  employed;  hence, 


(ISO) 


These  stresses  are  in  the  direction  of  the  laminate  reference  axes  X,  Y. 
The  stresses  in  each  ply  with  respect  to  the  natural  axes  of  the  ply  are 
obtained  from 


For  the  case  when  only  a  uniform  temperature  ciiange  is  present,  we  find  from 
equations  149  and  151  that  the  i^^-  ply  stresses  become 


(152) 


or  the  ply  stresses  in  its  natural  coordinate  systan  are  proportional  to  the 
differences  between  the  laminate  coefficients  of  thermal  expansion,  as  defined 
by  equation  147,  and  the  ply  coefficients  of  thermal  expansion. 


,A  ccm5)uter  program,  AC- 40,  has  been  avritten  to  perform  these  calculations 
for  4x4  for  any  crossplied  laminates.  The  program  listing  and  sample  data 
and  output  sheets  are  presented  in  another  section  of  volume  III.  Correlations 
of  this  theory  with  available  test  data  and  typical  design  curves  are  de¬ 
scribed  in  volune  I  of  this  report. 


SECTION  VIII 
JOINT  ANALYSIS 


SINGLE- LAP  ADHESIVE  JOINTS 


nie  derivation  of  the  single-lap  shear  joint  stress  distribution  equation 
for  relatively  flexible  adliesive  layers  is  shoi>m  in  reference  8  and  will  not  bo 
repeated  here. 


The  folloiving  assun5)tions  i\rere  used: 

1.  The  laminate  and  adhesive  are  isotropic  and  elastic. 

2.  The  laminate  deformation  is  cylindrical. 

3.  The  adhesive  is  relatively  flexible. 

4.  Hie  transverse  normal  strain  and  shear  strains  in  the  laminate  are 
negligible  in  comparison  iv’itli  tliose  of  the  adliesive. 

Tlie  results  are 


f" 

a 


where 


N 

> 

4L 


cosh- 


2^x 


^(1  +  300 


sinlr  0 


+  3(1  -ot) 


2(1  -  V-  )G  t 

xy  s 


E  t 
X  a 


I/: 


(153) 


a  = 


1  +  2y/r  tanh 


L  /K 
a/ 

'>  v 


3(1  -  V-  )N 
xy  X 


2E  t 

X 


Figure  45.  Typical  Single  Lap  Joint 


The  axial  stress  distribution  is  also  given  in  the  reference  but  will  not 

be  presented  here  since  it  was  not  used  in  our  analysis.  The  maximum  shear 

stress  occurs  at  x  =  t  L  /2: 

a 


fj  =  \  [^(l+3o0coth)?+  3(1-00]  ^ 

“a 


(154) 
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Substituting  equations  156a  and  156b  into  equation  155  yields 


•r*  =  cr 
?  X 

y  sin‘'0  +  -Hil- 
a  E 

/  2  -2 
(cos  0  -  v^y 

X 

- 

The  shearing  stress  is 

sin  </>  cos  (f> 

Using  Tresca's  maximum  shear  criteria,  we  find  that 


-  2 , 
sin  0 


(‘ 


DOUBLE-LAP  ADHESIVi:  JOINT 


1/2 


,,  Ea 

+  *^xy  g- 


\  Ea  2  ' 


2  2 

+  sin  0COS  0 


The  derivation  of  the  double-lap  shear  joint  stress  distribution  is  pre 
sented.  The  assun?)tions  are: 

1.  The  adherend  and  adhesive  are  isotropic  and  elastic. 

2.  The  adherend  displacements,  u.  and  u  ,  are  small  and  a  function 
of  X  only. 

.3.  Only  shearing  deformation  in  the  adhesive  is  considered. 


Figure  47,  Typical  Double  Lap  Joint 


Figure  48.  Typical  Load  Transfer  Through  Adhesive 
Considering  assumptions  2  and  3, 


and  the  constitutive  equation 

?  =  =  t;  (^2  - 


(158) 


(159) 


where 

G_  =  shear  modulus  of  the  adhesive. 

a 

Equilibrium  for  the  upper  element  requires 


The  incremental  elongation  of  the  adherends  can  be  expressed  as 


(160) 
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Cl61a) 


(161b) 


where  and  E2  are  the  Young's  moduli  o£  extension  of  the  adherends.  Noting 
assumption  2,  we  see  that  total  derivatives  can  be  used.  Differentiating 
equation  160  and  using  equations  161a  and  161b,  we  obtain 


^2^2 


0 


Equilibrium  for  the  joint  requires  that  the  condition  P  ®  +  P2  be 

satisfied;  hence, 


2G, 


(162) 


The  procedure  followed  is  that  of  reference  9.  However,  there  is  a  difference 
in  the  signs  of  the  last  two  terms.  The  solution  of  this  differential  equation 
is 


P^  =  A  sinh  wx  +  B  cosh  o/x  +  /3P  (163) 


1^ 

/3 


“Vi/ 


A,  B  =  undetermined  coefficients. 
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Tlie  boundary-  conditions  are 


(0)  =  P 


(164a) 


”1  C-a) 


=  0 


(164b) 


Solving  equation  163  for  the  undetermined  coefficients  yields 


(1  -^)  coshgjLa  sin]ia)x  +  (1  -jd)  sinh  wx  +yS 
sinli  ojL 


From  equation  (160) 


s  _  dPi 
^a  dx 


cosh  wx  -  (1  -/3)sinha>x  uP  (165) 

sinh  ojL^ 

Cl 


Practical  considerations  show  that 


ojLq  >  5 


0<)3<  1 

Thus  we  see  that 

(l-y3)  cosh  Lg 
cosh  wLg  s  sinh 
and  equation  (165)  becomes 


-s 


r  coshojL, 


^a  =  [iiHhZ; 

This  equation  attains  a  maximum  at  either  x  =  0  or  x  =  L^,  depending  on  j3. 
Evaluating  at  x  =  0,  we  get 


ft  (0)  =  (l-j8)a)P 

Cl 
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Evaluating  equation  (165)  at  x  =  L 


f!  (La)  =  (1-0) 


cosh^  <*^L„-sinh2  wL, 


sinh  wL, 


+  p  coth  uL  U)V 


'  iishr  ♦/’cothuLjjwP 


Therefore, 


fj  (I-a)  ~ 


Letting  P  =  Njj/2,  the  equations  become 


f|(0)  .  (1  oiLa  •  ^  -  (1  -fi  u,iY2 


PCl-a)  =  ^0,13  ^  =  ^a.N^/2 

3 


The  maximum  shear  stress  becomes 


a,  max 


=  IT  = 


where  )?*  is  the  larger  of  or  1-/? 


(166a) 


fl66b) 


(167) 


I 

I 


itelj.  ■.Wiiiw.' "  -  . . .  " 


DOUBLE-LAP  VARIABLE  STIFFNESS  ADHESIVE  JOIl^fT 

Hie  derivation  of  the  double-lap  two-adhesive  joint,  as  shown  by  fig¬ 
ure  49,  is  presented  in  the  following  paragraphs. 


P  - 


CT-tJ 


Figure  49.  Typical  Variable  Stiffness  Adhesive  Double  Lap  Joint 


The  procedure  used  here  will  essentially  follow  Szepe's  (reference  9).  How¬ 
ever,  his  equations  were  modified  in  reference  10  and  developed  in  the  previous 
section;  thus,  frcm  equation  162  and  figure  58  we  have 


l,xx  t  I  t,E, 
*  a  \  1  1 


1  .  2 


^2^2; 


=  -  2 


^a^2^2 


(168) 


The  solution  to  this  equation  is 


=  A  sinh  wx  +  B  cosli  wx 


(169) 


where 


w  = 


0  = 


'I  (vi  *  Vz) 


2  t^E^ 


2  ^2^2 
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Considering  the  three  regions  designated  as  one,  too,  and  three,  respectively. 


we  get 


=  A,  sinhoj^x  +  coshoj.x  +  /3P  0<  x<L  , 

11  1  1  1  1  al 


P„  .=  A-  sinh  W-X  +  B-  cosh  a)_x  +  flP  L  <x<L  . 
12  2  2  2  2  a^^  a2 


(170) 


P,_  =  A_  sinh  oj^x  +  B_  cosh  o)  x  +  flP  L  x  <L 
13  3  13  1  a2  £ 


noting  that  the  shearing  stress  in  the  adhesive  is  given  by 


=  P. 

a  1,  X 


fSi  =  cj  A  cosh  w,x  +  ct)  B  sinli  w.x 
a  11  111  1 


=  ^2^2  ^2^  ^  ^2^2  ^^2^ 


=  w,A_  cosh  w,x  +  w,B_  sinh  w^x 
1  :>  1  13  1 


(171) 


The  unknown  coefficients  A]^,  A2,  A3,  Bj,  83,  and  B3  can  be  obtained  from  the 
kno\vn  boundary  conditions. 


Pj^CO)  =  P 


’12  (^al) 


(172) 


**12  (*'a2)  “  '*13  (‘"02) 

°a2 


0 

Substitution  yields  the  following  equations 
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The  solutions  of  these  equations  are 
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where 


‘11  * 

sinh  w  L 

<^11 
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sinh  ca,L 
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S-,  «  sinh  <dLL 
21  2 


C-,  «  cosh  <j^l 

mX  to 


»  sinh  w_L 
Z  3.« 


C,-  »  cosh  6>-L 

22  2  a. 


S,,  »  sinh  W-L 
13  la 


C,  _  =  cosh  w,  L 
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11  21  ^2  21 
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X  *  c  c 

'^22  ^22  ^  13 
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^  ’  ^11  (^22  ^33 


0)2  21  '^23^ 


■  ^12  (^21  ^33 


wi  \ 

■  ;7^^22^23j 


The  loading  of  the  landnate,  P^^Cx)  and  P2W»  where  P^^W  ♦  ^2^^^  “ 
s 

the  adhesive  f  (x)  can  now  be  found.  To  check  these  equations,  we  examine  the 

dL 

case  idiere  a  single  adhesive  is  considered,  i.e. , 


P;^j(x)  for  Oix<L^j  ;  P13W  L^<  x  <L^ 


*■32  *  ■-al  ■  V 
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For  this  case,  the  equations  reduce  to  those  of  reference  3.  We  examine  next 
the  equations  for 


for  OsxsL^  and  =  0  ;  with 


For  this  case,  the  equations  aiso  reduce  to  those  of  reference  3. 


To  avoid  computational  difficulties  which  occur,  the  equations  can  be 
written  in  the  following  form.  Let  w=  jj—  . 
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Noting  the  following  relationships. 


(178) 


2  2 
cosh  X  -  sinh 

sinh  (x±y)  = 

cosh  (x±y)  = 


X  =  1 

sinli  x-cosh  y  ♦  cosh  x-sinli  y 
cosh  x-cosh  y  ±  sinh  x-sinJi  y 
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we  obtain  the  following: 


where 


Si  - 

■‘’41 

-  (l-SiCjjSjj 
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■^42  ^ 
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(179b) 


■*■  (^“^3  ^S-.l  I A 


11  55  ^2  66^ 
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11  12  24  -'3 
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-iS 

(179c) 


=  cosh  [(W2-w^)C^l 
C^2  =  cosh  ((W2-W^)C21 

=  cosh  [(<^2-^3  (C2  -q)] 
Cgg  =  coshtw^C^  -  (W2'^i3^23 

^6  '  («2-<^l)Cil 


=  sinh  ((W2-‘‘^_^3C^] 


S42  =  sinh  [(w2-W^3C2l 


S44  =  sinh  [(W2‘‘*i3(C2-C^3l 


S24  =  sinh  1^2(02-0^31 


=  sinh[W2Ci  '  <^^2'‘^1^^2^ 


It  should  be  noted  that  the  equation  leading  to  the  solution  of  A^,  B3  is 
still  almost  singular,  and  double-precision  arithmetic  must  be  used  on  a 
computer. 
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These  equations  can  now  be  used  to  determine  the  loading  in  the  laminate 
and  the  shear  stress  in  the  adhesive.  Figure  50  shows  a  typical  distribution 

of  The  discontinuities  of  the  shearing  stress  at  x/L^  =  0.25,  0.75  are 

violations  of  equilibrium  principles.  However,  considering  the  following: 

1.  only  six  free  constants  are  available  to  satisfy  boundary  conditions, 

2.  finite  element  solution  of  a  similar  single- lap  variable  joint  using 
an  energy  approach  (reference  12)  also  shows  an  abrupt  change  in 
stress, 

3.  the  thinness  of  the  adhesive  with  respect  to  the  joint  members, 

it  is  obvious  that  this  solution  is  basically  correct  although  the  actual 
stresses  would  be  more  like  the  dotted  line  shoim. 


p  =  Nx 
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Figure  50.  Double-Lap  Variable  Adhesive  -  Typical  Distribution  of  f®  versus  x/L, 
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APPENDIX 


FOlUiWORD 


Tins  appendix  contains  a  nimiier  of  computer  programs  tliat  were  generated 
during  the  theoretical  development  efforts  of  this  contract.  They  are  quite 
general  in  applicability  and  are  presented  here  separately  from  the  main  body 
of  this  volume  to  facilitate  their  usage.  For  eadi  program,  there  are  pre¬ 
sented  an  abstract,  the  input  data  format,  the  output  data  format,  and  an 
exantjle  problem,  'flie  following  is  a  list  of  tlie  program  nimbers  and  titles. 

No. 

Title 

AC- 2 

Extensional  and  Flexural  Constants  for  Balanced 

Symmetrical  Laminates  of  the  Form  [0,90/- 6/0 ,90]^ 

AC- 5 

H/C  Sandwicli  Panel  Stability  Under  In-Plane  Biaxial 

Loading 

AC- 7 

Analysis  of  Simply  Supported  Orthotropic  Sandwich  Flat 
Panels  Under  Pressure  or  Pressure  Plus  In-Plane  Loads 

AC-11 

H/C  Sandwidi  Panel  Stability  for  In-Plane  Shear  Loading 

AC- 12 

Hat-  and  Z-Stiffened  Panel  Stability  and  Stiffness 

AC-40 

N-Ply  Laminate  Theimoelastic  Analysis 
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■CGflPUrER  PROGRAM  AC- 2 


•  rEXIENSIONAIj.  AND-FLE^JRAL  .CONSTANTS  FOR: 
:BAM&fcSYMNiEfRIG^  iMHA-T^  OF  THE  KRM 

-^■■;(0,90/±^/0i90!^g 


ADVANCED  FILAMENTARY  COMPOSITE  MECHANICAL  PROPERTIES 
COMPUTER  PROGRAM  DECK  NO.  AC- 2 


ABSTRACT 


Fortran  program  deck  No.  AC-2  is  a  program  to  calculate  the  laminate 
flexural  and  extensional  constants;  the  honeycomb  (H/C)  sandi-/ich  parameters 
a.  and  and  the  magnitudes  of  the  flexural  rigidity  elements  D^.  and  D^^ 
relative  to  the  D-^  and  rigidities  by  using,  as  input  data,'^the  average 
mechanical  property  values  for  an  individual  unidirectional  lamina  and  the 
geometric  properties  of  the  total  laminate. 

The  program  is  somewhat  restricted  as  to  the  type  of  laminate  configura¬ 
tion  that  can  be  analyzed.  The  restrictions  are: 

1.  The  laminate  layup  is  symmetric  relative  to  its  midplane. 

2.  The  layups  must  be  balanced  with  respect  to  +B  and  -B  plies  (i.e., 
there  must  be  a  +0  ply  for  every  -B  ply) . 

3.  The  ply  orientations  only  include  one  angle  B  other  than  0  and  90 
degrees . 

4.  The  laminate  is  made  up  of  an  even  number  of  identical  sets  of 
ply  orientations. 

5.  Every  ply  is  assumed  to  have  the  same  nominal  ply  thickness. 

6.  Eveiy  ply  is  assumed  to  have  the  same  average  mechanical  properties. 

The  method  used  to  find  the  "exact'*  flexural  constants  differs  from  the 
current  state-of-the-art  method  in  that  no  summations  are  involved.  The 
summations  are  replaced  by  algebraic  total  sums.  Hiis  means  that,  for  a 
laminate  with  many  plies,  the  confuting  time  will  be  shortened  considerably. 
Also,  the  con5)uter  results  for  a  laminate  witli  many  plies  will  be  more  accu¬ 
rate  because  of  the  reduction  of  roundoff  errors. 

The  "exact"  flexural  constants  mentioned  come  from  the  standard  state-of- 
the-art  method.  The  program,  hcuvever,  also  selves  for  the  laminate  flexural 
constants  that  result  from  the  laminate  extensional  constants.  The  latter 
flexural  constants  are  bnly  approximations  and  consider  the  rigidities  D26  and 
D22  to  be  zero.  However,  if  a  large  number  of  plies  is  considered  (i.e.,  a 
large  number  of  sets) ,  the  "exact"  and  extensional  constant  methods  will  yield 
virtually  the  same  results. 
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PROGRAM  INPUT  DATA 


The  input  data  are  in  tv^ro  phases.  The  first  phase  (data  card  A)  consists 
of  giving  the  conputer  the  average  mechanical  properties  of  a  single  unidirec¬ 
tional  lamina.  These  properties  are  (longitudinal  Young's  modulus), 
(transverse  Young's  modulus),  (shear  modulus),  and  (Poisson's 
ratio) .  An  indicator  (A)  is  also  part  of  the  data  on  the  first  card.  Its 
purpose  is  to  tell  the  conputer  which  kind  of  average  mechanical  property  data 
are  being  used,  (i.e.,  A  =  1.0  for  average  conpression  data,  A  =  0.0  for 
average  tension-compression  data,  and  A  =  1.0  for  average  tension  data).  The 
correct  order  of  the  data  on  the  first  data  card  is  A,  E,  ,  E  ,  G^^  and 
The  input  format  is  (5F12.0).  An  exanple  phase  I  green  sheet  is  shown 
folloiving  the  text  in  this  section.  Note  that  all  five  input  values  need 
decimal  points. 

Ihe  second  phase  of  input  data  (remaining  data  cards  B)  consists  of  one 
data  card  for  each  laminate  configuration  (set  geometry)  to  be  analyzed. 

The  information  needed  on  each  card  is  IQ  (the  total  number  of  ply  sets) , 

H  (the  nominal  ply  thickness) ,  IPl  (the  number  of  inner  plys  in  each  set) , 

IPT  (the  number  of  ±0  plys  in  each  ply  set) ,  IP2  (tlie  nimiber  of  outer  plies 
in  each  set) ,  THETA  (the  ply  angle  associated  with  IPT) ,  IPSl  (indicator  for 
IPl) ,  and  IPS2  (an  indicator  for  IP2) .  The  format  for  phase  II  input  is 
(14,  F8.0,  314,  F12.0,  214).  This  means  that  the  integers  are  IQ,  IPl,  IP2, 
IPSl,  and  IPS2,  and  the  numbers  requiring  decimals  are  H  and  THETA.  An 
exanple  phase  II  green  sheet  is  shorn  in  this  section  following  the  text. 

Note  that  for  the  integer  numbers  (no  decimals) ,  the  format  calls  for  them  to 
to  be  against  the  right  side  of  the  allowed  width  indicated  on  the  sample 
green  sheets.  For  integer  numbers,  a  zero  is  placed  by  the  conputer  wherever 
there  is  a  blank  on  the  data  card. 

The  total  number  of  sets  indicator  IQ  must  always  be  an  even  integer 

(2,  4,  6  . )  because  of  the  assumed  symmetry  of  the  layup.  The  number  IPl 

is  the  total  number  of  0  or  90  degree  plies  nearest  the  centerline  of  symmetry 
of  the  laminates  as  shown  in  figure  1.  It  can  be  zero.  IPT  is  the  number  of 
±0  plies  in  each  set.  For  exanple,  if  IPT  =  3,  then  there  are  three  ±0  pairs 
in  each  ply  set  (i.e.,  six  plies).  It  can  be  zero.  IP2  is  like  IPl  except 
that  it  is  the  total  nunber  of  0  or  90  degree  plies  in  each  set  tliat  are 
to^^?ard  the  outside  or  surface  of  the  laminate.  It  also  can  be  zero.  THETA 
is  the  plus  and  minus  angle  of  the  middle  plies.  The  program  assumes  that 
for  an  input  of  +0,  the  +0  ply  starts  closest  to  the  centerline.  If  the 
opposite  is  desired,  then  the  THETA  value  must  be  input  as  a  minus  THETA. 

IPSl  indicates  the  angle  of  the  inner  plies.  If  it  is  zero,  then  the  inner 
plies  of  each  set  are  0  degrees,  and  if  it  is  90,  then  the  inner  plies  are  at 
90  degrees.  It  cannot  take  any  value  other  than  0  or  90.  IPS2  is  like  IPSl 
except  that  is  refers  to  the  outer  plies  in  each  set. 
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0*  OR  90* 


MIDDLE  PLIES 

2pgt* 

I _ 


0*  OR  90* 


INNER  PLIES 


ONE  SET 


OUTER  PLIES 


t*  -  NOMINAL  PLY  THICKNESS  (H) 
tg  =  SET  THICKNESS  =>  t*  (  l  +  2p^  +  2) 
tt  =  total  laminate  THICKNESS  -  qtg 
q  «  NUMBER  OF  SETS  (iQ) 


Figure  1.  AC2  Laminate  Mechanical  Properties  Program  Figure 
for  Ply  Set  Orientation 


o  o  o  r»  on  o*r>  oooo  nonoooono  non  ooooo<^onooooooo 


.PRflRRAM.ACZ  EXTENSIONAL.  AMO  FLEXURAL  CONSTANTS_j:flii. BALANCED i 
SYMMETRICAL  LAMINATES  OF  THE  FORM  (0/+-THETA/90 )N5 


THE  INPUT  OATA  IS  AS  FllLLIIHS . 

THF  FIRST  DATA  CARD  HOST  CONTAIN  THE  LAMINATE  ELASTIC  CONSTANTS  IN 

_ THE  FOLLOUING  FORMAT.. .(bFl^.O)  .  _  _  _ _ 

THE  LAMINATE  ELASTIC  CONSTANTS  ARE  A , EL, ET , GLT , UL T  WHERE . 

A  IS  A  CONTROL  CARO  THAT  INDICATES  WHICH  TYPE  OF  ELASTIC 
CONSTANTS  ARE  TO  BE  USED.  A=-1.0  IS  FOR  COMPRESSION 
LOADINGS,  A=0.0  IS  FOR  TFNS ION-COMPRESS  ION  LOADINGS,  AMU 
A=1.0  IS  FOR  tension  LOAOINGS. 

_  EL  IS  THE  LAMINA  LONGITUDINAL  YOUNGS  MODULUS 

ET  IS  THE  LAMINA  TRANSVERSE  YOUNGS  MODULUS 
GLT  IS  THE  LAMINA  SHEAR  MOOULUS 
ULT  IS  THE  LAMINA  POISSONS  RATIO 

THE  SECOND  AND  FOLLOWING  DATA  CAROS  MUST  CONTAIN  THE  SET 

_ DESCRIPTIONS  IN  THE  FOLLOWING  FORMAT. . . { I A , FB .0, 314 , F12 .0, 214 )  _ 

THE  SET  OESCRIPTIONS  INCLUDF  0,H, i P 1 , I PT , I P2 , THETA, I  PS  1 ,  AND  IPS? 

WHERE . 

10  IS  THF  TOTAL  NllURFR  OF  SETS  IN  THE  LAMINATE  (  IT  MUST  ALWAYS 
BF  AN  EVFN  INTEGER) 

H  IS  THF  NOMINAL  PLY  THICKNESS  (FLOAT.  PT.I 

_ IPl  IS  the  number  OF  INNER  PLYS  IN  EA(.H  ID  SET  IN  THE  7.EKU  OR  . 

OEGRFE  DIRECTION  DEPENDING  ON  IPSl  (INTEGER) 

IPT  IS  the  NUMBEIi  OF  MIDDLE  PLIES  IN  EACH  O  SET  IN  THE  PLUS  AND 
MINUS  THFTA  OEGRFE  DIRECTION  (INTEGER) 

IP2  IS  THF  NUMBER  OF  OUTER  PLIES  IN  EACH  ID  SET  IN  THF  ZERO  OR 
degree  direction  depending  on  IPS2  (INTEGER) 

_ THETA_IS  THF  PLUS  AMD  MINUS  ANGLE  OF  THE  P2  PLIES  .(FLOAT,  PT.)-.  ..  , 

IPSl  IS  AN  INDICATOR  TO  DETERMINE  IF  THE  ANGLES  ASSOCIATED  WITH 
PI  ARE  ZERO  OR  NINETY  DEGREES.  IF  IPS!  IS  ZERO  THEN  THESE  ANGLES 
ARE  ZERO  DEGREES.  IE  IPSl  IS  90  THEN  THESE  ANGLES  ARE  90  DEGREES 
IPS2  IS  THE  same  AS  IPSl  EXCEPT  THAT  1PS2  REFERS  TO  P2 


DIMENSION  0(A),0PR(4,3I,UC(3) 

READ(5,  9991  A, EL, ET, GLT, ULT 
999  FORMATISFIZ.O) 

2000  READ(5,1000)  IO,H, I P 1 , I PT , I P2,THETA, I PSl , I PS2 
1000  FORMAT! 14, F8. 0,314, F12.0,2I4) 


_ IF(A)7,5,3 _ _ _ _ _  _  _ _ 

3  HR1TE(6,4) 

4  FORMATCI... LAMINATE  ELASTIC  CONSTANTS  FOR  TENSION  LOADINGS') 

GO  TO  14 

5  HRITE(6,6) 

6  EDROATI • 1... LAMINATE  ELASTIC  CONSTANTS  FOR  TEN-COMP  LOADINGS') 

GO  TO  14  _ 

7  WRITF('6,*n) 

8  FORMATCI. ..laminate  ELASTIC  CONSTANTS  FOR  COMPRESSION  LOADINGS') 
14  WR1TE(6,9)  EL‘,ET,GLT,ULT 

9  FORMAT! >0  EL='E12.5,5X«ET='E12.5,5X'GLT='E12.5,5X'MULT='F8.4/// ) 


C 

C _ 


0  =  10 
P1=IP1 
PT=IPT 
P?=1P2 
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[1  Hlvtltt" '''tilt '"!(>■  p'i'r  ,  I  f  III' P'l'I 


RK=PI+2*PT+P2 
_ XM»n«RK 

■ri=xn/2-rk  ■ 

R2=R1+P1 

R3=XN/2-P2 

T=XN>»H 

IPT2=2*PT 

C 

1=0 

PM»P1 

R=R1 

100  1=1+1 

V=n*PN/8={ 2* tPM+1 )*(2*PN+6*R-3*XN-?)+3*{2’»R-XN- 

_ 1PM+8<=R-3»XN+RK=>(  1+01  1  +  1  1 

IFtI-21  110,120,130 
110  DC ( I  1 =V 
PM=2«PT 
R=R2 

GO  TO  100 

- 1Z0J)CJXI=V _ 

PN-P2 

60*^70  100 
130  0C( 1 1=V 
C 

C=C.nsn<THeTAl 
S=S1MD( THETA  1 

_ C?S?=(C2Si«*2 _ 

S2=S=»*2 

C2=C*’»2 

C2TS2=C2*S2 

C4=C<=*A 

54=5*44 

_ Cf  flS4sC4±SA _ 

C 

UTL  =  (6T/ELI*llLT 
PS1=1-ULT>»UTL 
011=EL/PS1 
012=011*UTL 

_ Q22=ET/PSJ _ _ _ _ 

066=GLT 

C 

OPR ( 1 ,2 )=P1 i*C*+2*(gl2+2»066 ) *C2S2 +022*54 
OPRl?, 11=012 

OPR (2,21= (01 1+022-4*066 1*C2S2+012*(C4+S4) 

_ pPR(2, 31=012 _ • _ 

C'PRl 3, 2 1=01 1*54+2* (012+2*066 )*C2 52 +022*C4 
0PR(4, 11=066 

OPR (4,21= (01 1+022-2*012-2*066 l*C252+Q66*(C4+S4l 
0PR(4, 31=066 

1F(1PS1-901  140,145,145 

_ litfl..OPR.{ij.Il=QlJ _ I - 

•0PR(3, 11=022 
GO  TO  150 
145  QPR(l, 11=022 
0PR(3, 11=011 

150  1F(IPS2-901  155,160,160 

_ 155_QPR  (1,3  }.=Qil _ 

0PR(3, 31=022 
GO  TO  165 
160  OPR (1,3 1=022 
0PR(3, 31=011 


l**2+RK*(0+21*(3>i‘ 


4 
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165  CONTINUE 


AN=.6666666*(T/XN)*>t>3.0  ^ 

on  170  I»l,4 
0(I)=0.0 
on  170  J=l,3 

170  0»  n=0<  I  )+OPR(  1,J  )*DC(  J)*AN 

’*'‘'**3*5+ ( 012-022+2*066  )=>C*S««3 

- ^2*T^l«ll^Pl?~2*0f^)*C*.S**3  +  (O12r.^^  _ 

D{6)=026T*CnN 

- RAT1QL=D15)/D(11  . 

RAT102=D(5)/n{2)  ■  '  - -  '  — . - 

RATtn3=0{61/n( 1 ( 

RATlnA=DI6)/n(2) 

A1  =  P1«:0PR<  1,  1 )  +2*PT*0PR(  1  ,2)  +P2vfipsM  3, 

- "*^=^1*^'’'*<2.11_+2*PT*  QPR(2,2)  +JP2*QPR(2,a, 

A2  =  P1’»0PP,(3,1)  +2*PT*0PH(3,2)  +P2*npft  ( 3,  3  ,  '  "  . 

A6=P1*0PH(A,1 )  +2*PT*0PR(4,2)  +0PR(4,3)*P2 

11YX=A12/A1 
UXY=M2/A2 
_ f'XY=A6/RK 

6X=( Als<A2-Al2«*2') /A2/RK  ’  ’  •• 

EY={A1*A2-A12**2)/A1/RK 

AI.F=(EX/FY)*«.5 

PS  IfstEX’JFY)**. 5/(1  .0-UXY*!IYX) 

B  =  1-I)XY*UYX 

- AI  l=FX«(  1.0+r-Y/ex+2,0*IIVXI/B 

A12sGXY-Fx*lJYX/B  ' 

KRITF{6,575) 

FORMAT! 'O,,, SET  OFSCK IPT ION • 1 
WRITE(6,576) IO,H,T 

OF  PLIES  WITH  A  NOMINAL  PLY 

P^INrHFS  V)  lahinatf  thickness  IS'F 


575 

576 


THICK. 

10.4' 


577 

578 


579 

580 

585 

590 


C 

C 


600 

610 


IF(IPS1-90J  577,579,579 
..WRITF(6t5JfiJ.IPl . . . . 

FORNATC  NUMRER  OF  INNER  PL  IFS*  A~T  0  DECREES  >•  1 3)  ”” 

GO  TO  585 
HRITE(6,5R0J  IP! 

FORHATC  NUMBER  OF  INNER  PLIES  AT  90  DEGREES  ='I3> 
WRITF(6,590)  theta, IPT2 

F.;rmat(»  ...number  of  middle  PLIES. at  (+  and  -»*F7.2*  degrees 

IF(IPS2-90)  600,620,620 

HRITE(6,610)IP? 

Gn'^Tn^635  degrees  =*I3//I 


.=2131. _ 


620  WRlTE(6,630)r,P2 
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f  ^  630  F0RI*6T(*  NUMBER  CF  CUTER  FLIES  RT  90  DEGREES  -MS//! 

f  635  WRITE«6,66CJ_ 

5  640  F0BHRT<*0...L*HIN*1E'  FLEXURAL  CGNSTANTS’/l 

I  C 

f:  URlTE<6t650)  0 1 1 ) t0(2) t0( 31 1 A1 . A12 10(4) tOIS) *0(6) *62 1 AE 

650  F0RHAT(*011  **E12.5*3X<0I2  « •E12.5,3X*C22  »*E12.5,7X*A11  -•E12.5t 
E  15.'<*A12  «*El2.5/‘066  >«E12.5f3X*D16  -•Et2.5*3X*C26  «'E12.5*7X'A22  « 

;  _  2«EI2.5.3X«A66  <'E12.S//) 

.  {.  - 

I  tiRITE(6,660) 

r  660  FORHATCO... RATIOS  BETWEEN  FLEXURAL  CONSTANTS*/) 

i  C 

I  WRlTG(6f670)  RATIOl tRAT 102 tRAT 1C3,RAT 104 

?  ^0  F0RHAT(*016/C11  3*E12.5*SX*016/012  >*E12.5t5X*C26/011  >*E12.Sf5X*0 

fc.  126/012  **E12.5//)  ■  ■■ 

6  C 

£  WR|TE(6,6eC) 

S  660  FORNAT(*C...LAFINATE  EXTENSIONAL  CONSTANTS*) 

f  C 

f  _  )|R1T6(6,13)  EX*EytGXY*UXy,UVX.ALF.PSlI»AIl*At2 

S  13  F0RHAT(*0*4X*EX»*E12.5,5X*EY=*E12.5,5X*CX¥a*El2.5,5X*NUXy»*F7V4,5X 

'  l*MtYX**F7.4/*  * IX* ALPHA-SORT (EX/Ey)>*E12.5. 13X*FS1-S0RT( EX«EY)/(1- 

-  2PUXV»HL'VX)=*E12.5/*  *3X*A11-EX*( UEY/EX42«HUyX)-*E12.5tSX*AI2-6XY- 

#  3EX*HUyx/(l-HUXy«MUYX)>*E12.S//) 

c 

_  BC-T«*3/12/B _ 

6iia-bc«ex  ■ 

012A-UyXaEX*BC 

r-  _ 022A-BC«EY _ 

066A*GXY«e«8C 

-  C 

WRITE(6*7CC) 

f:  7C0  FORMAT! *C...LAFtNATE  FLEXURAL  CONSTANTS  FROM  LAMINATE  EXTENSIONAL 

ICONSTANTS*/) 

c  . ' 

HRITE(6*71C)  CtlA*012A*022A,066A 

710  F0RMAT(*011E  **E12.S*5X*012E  >*E12.5*5X*D22k  -•E12.5*5X*066E  **E12 
;  1.5) 

I 

I  GO  TO  2000 

f  END 
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EXAMPLE  PROBLEM 


AC- 2  LAMINATE  FUEXURAL  AND  EXTENSIONAL  PROPERTIES 

Consider  the  three  laminates:  [0/±45/0]g,  [02/^45/902] g»  and  [0/±60]g. 
AC- 2  is  used  to  find  the  elastic  and  flexural  cmstants  from  the  single-ply 
properties  which,  for  room  tenperature,  are: 

=  29.9  X  10^  psi 

Ej.  =  2.71  X  10^  psi 

=  0.70  X  10^  psi 

'll  = 

The  data  sheets  for  these  three  laminates,  as  well  as  the  output  printouts, 
are  shcnvn  at  the  end  of  this  section. 


L«HINATE  ELASTIC  CONSTANTS  FOR  TEN-COHP  LCAOINGS 


LAMINATfi  ELASTIC  CONSTANTS  FOR  TEN-COMP  LOADINGS 


. . . . . . 
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p^l.  ;st^iLiTY  i^er  in-plane  biahal  loading 


r 
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ADVANCED  FILAMENTARY  Ca^iPOS^rE  MECHANICAL 
PROPERTIES  COMPUTER  PROGRAM  DECK  NO.  AC- 5 


ABSTRACT 


Fortran  program  deck  No.  AC- 5  is  a  program  to  calculate  the  general 
stability  of  a  honeycomb  (H/C)  panel  with  orthotropic  core  and  face  sheets 
under  biaxial  loading.  The  loading  ratio  0(0  =  Ny/Nx)  can  be  zero,  plus,  or 
minus.  Thus,  AC- 5  will  handle  panels  under  uniaxial  compression,  biaxial 
tension-compression,  and  biaxial  con|)ression-conpression.  A  listing  of  AC-5 
is  shown  following  the  text  in  this  section. 


PROGRAM  INPUT  DATA 


The  input  data  are  shown  on  the  sanqjle  data  sheets  in  this  section.  The 
data,  except  for  the  title  card,  are  read  in  with  a  6F12.0  format  per  card. 
The  definitions  of  the  input  data  areas  follow: 


Card  A: 


TIT 

Card  B: 
'INA0B 
ALF 
TNN 
TOM 
TOPHI 

TOV 


Any  desired  one- card  title  in  columns  1-72. 
If  no  title  is  desired,  use  a  blank  card. 


Total  number  of  (a/b)'s  to  be  considered;  maximum  of  20. 
Alpha  =yjE^/'E^  (fran  program  AC- 5). 

Total  number  of  half-waves  in  the  y- direction;  no  limit. 

Total  number  of  half-waves  in  the  x-direction;  no  limit. 

Total  number  of  load  ratios  0;  maximum  of  20 
(from  program  AC- 5). 

Total  number  of  Vs;  maximum  of  20. 


Card  C: 

All  =  Invariant  1  =  E^(l  .  E^/E^  *  2 
can  be  found  from  program  AC- 5. 

AI2  =  Invariant  2  =  G^  -  E^  *yx4^  '  *'xy'yx)”’  found 

from  program  AC-5. 
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r.  =  Grthotropic  core  shear  nxDdulus  ratio; 

TNPSI  =  Total  number  of  tp;  maximum  of  20. 


Card(s)  D: 


A0B(I)  =  Panel  aspect  ratios  (a/b)^.  There  are  six  of  these  to 
a  card  tip  to  a  maximum  of  20  four  cards). 


Cardfs)  E: 

PHI  (I)  = 


Load  ratios, 0=  Ny/N^^.  There  are  six  of  these  to  a 
card  15)  to  a  maximum  of  20.  These  cards  are  to  be 
left  out  if  TMPHI  =0.0. 


Card(s)  F: 


Core  shear  parameters,  \diere  V  =  (c  +  t^)/^2G'  b^j 

There  are  six  of  these  to  a  card  up  to^a  maximum  ^ 
of  20. 


Card(s)  G: 


PSICD  = 


Modulus  ratios  found  from  AC- 2  where  JB  E  /(l-v  v  j. 
There  are  six  of  these  to  a  card  up  to  a  ^  ^  ^ 

maximum  of  20. 


An  input  data 
simply  stacked  one 


flow  chart  is  included.  Note  that  additional  cases  are 
behind  the  other. 
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CARD  C 


CARD(S)  D 


READ  V,.  V^.  -  -  - 


READ  ■  ■  ■  ho 


CARD(S)  E 


CARD(S)  F 


CARD{S)  G 


-NEXT  CASE 


OUTPUT 


The  output,  as  shown  in  this  section,  gives  the  title,  the  core  shear 
ratio  R,  gamma,  beta,  alpha,  psi,  and  V.  It  gives  and  ky  for  uniaxial 
compression  in  the  x  and  y  directions,  respectively,  for  zero  phi  (phi  =  N 
and  any  other  values  of  phi  as  given  in  the  input  data.  ^ 


O  U  O  O  O  u'o  UOOOOOOOOU  OlO  u  o 
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PROGRAM  ACS  ..'i  H/C  PANEL  STABILITY  UNDER  IN-PLANE  BIAXIAL  LOADING 

THE  INPUT  DATA  AND  SELECTED  CflNSTAMTS  ARE  AS  FOLLOWS .  . 

All  =  INVARIANT  1  =  EX/(1-UMXY<=UMYX)*(1+EY/EX+2*UMYX) 

AI2  =  INVARIANT  2  =  GXY  -  UMYX=>6X/ 1  l-UMXY*UMYX ) 

_ ALF  =_SnRT(EX/EY)  _ _  _ _ 

beta"  =  H  =  ALF*UHYX  +  2'<'GAMMA 
Oil  =  FXwTF#C**2  /  12.*  (  l-UMXY=»imYX  ) ) 

022  =  nil*EY/FX 
066  =  TF*C**2*bXY/2 
OOX  =  C*GCX 

_ aQY._=..C*GCYL _ _ _ _ 

GAMMA  =  G  =  GXY«{ 1-UMXY«I1MYX)/S0RT(EX*EY1 
lA  =  TOTAL  NUMBER  OF  A/B 

IC  =  TOTAL  N'lMBcR  OF  HALF  WAVELENGTHS  IN  Y-OIRECTION 

ID  =  TOTAL  NUMBER  OF  HALF  WAVELENGTHS  IN  X-DIRECTION 

IE  =  TOTAL  NUMBER  UF  PHI,  WHERE  PHI  =  NY/NX 

_  IF  r  TOTAL  NUMBER  OF  V  WHERE  V  =  P I  *»2»S0R  T  ( D1  1  *022  I  /  ( B*<=?*I)UX  1 

PST  =  SnHT{EX*FY)/( 1-UMXY«UMYX) 

R  =  OOX/OOY 

DIMFNSIOM  AOR1201,  PHI(20I,  V(20)»  AK(20),  APSI(20),  TITI18), 
1AAK(20) ,AAKAY(20) ,  AKAY(20) 

1  REAP  >  (TITII  )  ,  I  =  1  _ 

100  FORMAT!  1  baa") 

RFAOlStllO)  TNAOB,  ALF,  TNM,  TNM,  TNPHI,  TNV,  All,  A12,  R,  TNPSl 
110  FOHMAT(6F12.0) 

A^AI  F 
IAsTNAUB 

_ IC  =  T_MN _  _ _ _ _ 

jOsTHM””  - 

HEAn(5,llO)  (AOBII) ,1=1, lAI 
IF=TNPHI 

IFIIEI  10H,10H,10S 
105  RFA0(5,110)  (PHK I ),I=1,IE1 

lOR  1P=TMV _ 

READ(5,110)  (V(n,I  =  l,IP) 

L2=TNPSI 

_ ^READlfttllO)  ..lAPSl(I),I=ltL2I  ....  .  _ 

DO  300  L3=1,L2 
PSI=APSI(L3) 

GAMMA=  All/(2.*psn  +  AI2/PSI  -.S*A  -.5/A 

HFTA=  l,5*An/PSI  +  2.*AI2/PSI  -1.5/A 

I  _ R  =  BFTA _ _ 

■  (;=GAHHA . 

i  HR1TF(6,120) 

120  FORMAT! '1  BUCKLING  COEFFICIENTS  FOR  SIMPLY  SUPPORTED  ORTHOTROP 

;  lie  SANDWICH  PLATES  -  FLAT* 1 

WRITE!6,1001  !TIT!I),I=l,lfl) 

_ HRlTF!6,nOI  R,  G,  B,  A,_PSI _ 

130  FORMAT! •0*SX'G"CX/G* •CY=*F5. 2,5X*GAMHA=* F8.A,5X 'BETAS 'Fa. 3, 5X* ALP 
1HA=«F6.3,5X*PS1='F11.0) 

DO  300  1=1,  IP 
WH1TF!6,135)  V(I) 

135  FORMAT! lAHO  DATA  FOR  VsFB.A) 

_ IfULE  l_l3.8.»JA£b  138 - 

13B  WR!TE!6,139)  IPHI ILL ) ,LL=1, I E ) 

139  FORMAT! •  '26X 'PHI*' F5. 1 ,6X' PHI = 'FS.l ,6X 'PHIs'FS. I »6X' PHI='F5.l .6X • 
1PHI»*F5.1)  _ 

14C  continue 
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142  FnRMAT( «0«2X • A^BJ 5X «KX' 6X»K Y«5X *KX« 1 >/KY< 1 1 «4X 'KX (2 ) /KV (2 » «4X 'KX 1 3 
1 »  /iCY  (  3  r*  4X  'KxTaI  /KY<  4 )  '^X  '“kX  (  5» /KY  <5 » •  ) 
l)n  300  J*1,IA 
APl=AnH( J) 

AB?=ABi*>'2 
AR3=AB1  !!*3 

_ AB4=ABi5*4 _ _ _ _ 

AAKX=  I.OK+IA 
AAKY=AAKX 
IF(IB)  145,150,145 
145  ntt  14B  J.l  =  l  ,IF 
AAKAV(  J.I)=aAKX 
14M  AAK(.1.I)-AAKX 
ISO  f.MN  I  INlll- 

0(1  250  M=l,in 
OM=H 

nM2=nH>s*2 

OM3=nM*c3 

nH4=nH7->g2 _ _ _ _ _ 

00  250  N=1,IC 
ON=N 

ON2=ON**2 

flN3=nN4(|M2 

ri(>i4=nN2*«2 

Alls  (A*nM3/AB3  +nM*(lN2=»B/ABl)*(-l. ) 

A12=  1.  ♦  0M2«A*V(  11/AH2  +nN2«r,*v<  1 ) 

A13  =  (-l.  n»K*(B-G»/ABl  )  _ _ _ _ 

A21  =  (nN3/A“+nfJ*0MS«2wH/AB2) 

A22=  A13/R 

A23=  l.+B  a((3N2/A  +(;«()M2/AB2)«  VCD 

A31  =  (-l.  )=!=(  A*nM4/AB4  +2.«B*OM2sON2/AB2  +  0N4/A) 

A32=  -OH  /ARl«  iriM2«A/Ab2  ♦0N2«B»*V(n 

_ A33=_nNwR'»V(  I  )<=(()M?/A_*'  ClH2£i»^B2.l _ 

AP31  =  (A'3“3‘'»('Air=!'A22-A12»A2U  +A32*< A13*A2 l-All*A23 )  I /  (A12#A23  -A22 
1«A13) 

AP=  AP31-A31 
AKX  =  Ai‘/(llM2/Ae21 
AKY-AP/0N2 

fFIAKX-AAXXI  160,170,170 
160  AAKXsAKX 
J70  CONTINUE 

iFfAKY^A/'KY)  I75iI80.160 

175  AAXYsAKY 

IHO  r.dUlIMUH 

IFdFl  1RS,2H0,185  _ 

Th?" (in  195  C=1  ,IF' 

AK(l.)=AP/(  nH2/AB?  +(IN2wPHl(LH 
IF(AK(L)-AAK<L))  190,195,195 
190  AAK(Lt=AK<L) 

195  CONTINue 

on  19H  L  =  1,IF  _  ..  - - - - - 

AKAY(L)=AK{Lr«PHl  (Cl 
IF(AKAY(L»-AAKAY(Ln  197,198,198 

197  AAKAY(L1=AKAY(L) 

198  CONTINUE 
250  CONTINUE 

280  WR1TE(6.290)  A0B(J1.  AAKX.  AAKY.  lAAK  IL  IjAAl^  AY  I L  >  ,L  =  l  1 1 E ) - 

290  FORMAT!*  •F6.2,2F8.2,  5(FR.2'/*F6.2  )) 

300  CONTINUE 
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EXAMPLE  PROBLEM 

AC- 5  H/C  PANEL  STABILITY  UNDER  IN-PLANE  BIAXIAL  LOADING 

Hie  three  laminates  to  be  used  in  the  exanple  are  [0/*45/0]g,  2S’ 

and  [0/-60]j,  at  room  temperature.  The  values  for  a.  All,  AI2,  and  E  ,  E  , 


G  ,  V  ,»  and  v  can  be  found  from  program  AC- 5. 
xy  xy  yx 

The  gecmetry  of  the  panel  is  given  as 


X  y' 


c  =  O.S  inch, 


a/b  =  1,  2,  3 


b  =  15 


R  =  D  /D  =2.0 
qx  qy 

^  =  0,  .1,  .4,  1.0,  n,  m  =  1,  2,  ---6 

2.0,  10.0 


The  material  properties  of  the  core  are  such  that 
G'  .  =  126,000  psi  (8.1  Ib/ft^  core) 

Now,  for  [0/f45/0]g  laminate,  t^  =  8  x  .0052  =  0.0416  inch 

Qf=  1.815,  All  =  33.884  x  10^  AI2  =  0.12862  x  10^,  ^'=  10.794  x  10*^ 

„  _2.  /  ,2  (9.8690)  (10.794  X  10^)  (o. 0416)  (o. 500  +  0.0416) 

V  =  IT  ^tj,  (c  +  t^l/Zb  b  =  - ‘ - ‘ - rr-7 — - 

“  ’(o.I26xlO-'')(ls)^ 


=  0.04233 


For  j02/*'*5/902l2S  t^  =  24  x  0.0052  x  0.1248  inch 

a  =  1.00,  All  =  33.884  x  10^,  AI2  =  0.12862  x  10^,  4' =  13.971  x  10^ 

Y  =  (9.86961  (13.971)  (0.1248)  (0.6248)  ^  g 
2(0.126) (15)^ 


For  [0/-60]g  laminate. 


t^  =  6  X  0.0052  =  0.0312  inch 


a  =  1.00,  All  =  33.884  x  10^,  AI2  =  0.12862  x  10^,  4-  =  12.771  x  10^ 

y  ,  (9.8696) (12.771) (0.0512) (0.5512)  ^  g  53^84 
2(0. 126) (15)^ 
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COMPUTER  PROGRAM  -AG-7 

ANALYSIS  OF  SIMPLY  SUPPORTED  .  ORTliOTROPIC  SANDWICH  Fi.Af 
PANELS  -UNiJER  PRESSURE  OR  PRESSURE  PLUS  IN^PLANE  LOADS 


AC-7 


ADVANCED  FIL^MIiNWY  COMPOSITE  MECilANICAL  PROPERTIES 
COMPUTER  PROGRAM  DECK  NO.  AC -7 


ABSTRACT 


Tliis  program  calculates  the  deflection,  moments,  and  sliears  in  a  flat 
orthotropic  honeycomb  (H/C)  panel  under  pressure  plus  in-jiL'ine  loads,  given 
material  and  geometry  parameters.  Tlie  metliod  is  a  series  summation,  and  the 
series  terms  as  well  as  the  series  sum  are  i>rinted  out  to  enable  tiie  user  to 
monitor  convergence. 

The  in-plane  loads  can  be  compression-cofiipression,  tension-comi^ression,  or 
tens ion- tens ion.  For  the  first  two  cases,  part  of  the  imput  data,  tlie  buckling 
coefficients,  can  be  found  from  program  AC -5.  For  the  latter  case,  tension - 
tension  buckling  coefficients  are  found  as  noted  later  in  tins  program 
description. 


PROGRAM  INPUT  D.ATA 

The  program  uses  as  input  che  results  of  programs  AC -2  and  AC- 5,  as  well 
as  other  items.  Tlie  input  items  are: 

Card  A: 


TIT  =  Any  desired  one-card  title.  The  title  is  in  columns  l-“2. 
If  no  title  is  desired,  a  blank  card  must  be  used. 


Card  B: 


'DiA0B  =  Total  number  of  (a/b)'s  to  be  read  in  (ma,xLTium  of  20). 

.•\LF  =  Tlae  alpha  volume  from  program  AC -2  (or  =  -JUx/Ey) . 

TNM  =  Total  number  of  half-waves  in  the  X-direction  (no  limit) . 

T.W  =  Total  number  of  core  shear  influence  parameters  V  to  be 
read  in. 

•All  "  Invariant  1  from  program  AC-2.  It  is  E^(l  + 

2v  1//1  -  V  v  1. 
yx//\  xy  yx/ 

AI2  =  Invariant  2  from  program  AC-2.  It  is  G^  -  E^i/^^/^l  - 

V  V  \* 

xy  yx 
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Card  G: 


R  =  Ratio  of  N  to  critical  N  /N’ ./N  plus  for 

X  X  X  \  X  X^  C  i*' 

compression  and  minus  fur  tension. 


R  =  Ratio  of  N  to  critical  N  {N  />i  plus  for 

y  y  y  V  y  y,  cr; 

com^iression  and  minus  for  tension. 

R  =  Core  shear  modulus  ratio  G’  /G’  . 

cx  cy 

TNPSI  =  Total  number  of  ip  to  be  considered. 

TNX0A  =  Total  number  of  x/a  to  be  considered. 


TMY0B  =  Total  number  of  y/b  to  be  considered. 


CardCs)  D; 


.40B(13  =  Panel  aspect  ratio  (a/b)^. 


.'lAKXCl)  =  Buckling  coefficient  K  from: 

X 


1.  .‘\C*5  corresponding  to  given  value  of  0  if  is 
compression 


2.  K  =  N  ( D~~  if  .\  i 
X  x\7r//  11  22  X 


IS  tension. 


A-VKY(l)  =  Buckling  coefficient  from: 


1.  -AC -5  corresponding  to  given  value  of  0  if  is 
compression.  ^ 


2.  K  =  .\  D.,.,  if  N .  i 


y\7r 


11  22 


IS  tension. 


Tliere  are  two  sets  of  these  parameters  to  a  card  (maximum  of  20 


sets) . 


Card(s)  E: 


V(I)  =  Core  shear  parameters  where  V  =  )r“^t^^c  + 

TJiere  are  six  of  these  to  a  card  (maximum  of  20) . 
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Card(S)  F 


APSI(I)  =  Modulus  ratios  4>  found  from  AC -2,  where  $  =  - 

1/  w  ).  There  are  six  of  these  to  a  card  (maximum  of  20). 
xy  yx/ 

Card(S)  G; 

X0A(I)  =  Ratio  or  point  on  X-axis  where  results  are  to  be  given  to 
the  lengtJi  of  the  panel,  a,  in  the  x-direction.  There  are 
six  of  these  to  a  card  (maximum  of  20) . 


C:ard(s)  H: 


Y0B(I)  =  Ratio  of  point  on  Y-axis  where  results  are  tc  be  given  to 
the  width  of  the  panel,  b,  iii  the  y -direct ion.  There  are 
six  of  these  to  a  card  (maximum  of  20) . 


Sample  input  data  sheets  are  sliown  in  this  section  folloicing  the 


text. 
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CARD  A 


1  TN  a/b,  O',  TNM 

TNV,  All,  AI2  - 

R  ,  R  ,  R,  TN 

X  y 

TN  x/a,  TN  y/b 

f~ - - - ' 

READ  NEW 
CASE 


CARD  B 


•CARD  C 


CARD(S)  D 

.  CARD(S)  E 

-CARD(S)  F 

-CARD(S)  G 

-.CARD(S)  H 


additional  cases  may  be  stacked  one  behind  tlie  othei 


OlJTPUl 


Iha  output,  Shown  later  in  this  secti™  gives 


T  n’  r  r  a/b  V  r  a  The  series  terms  are  printed 
toT-  ;  ""I  3’  S^  -'-  TO  for  rte  values  of  x/a  and  y/b  given.  The  last 
Lrm  in  the  sU  Un  be  used  to  find  the  various  quantities,  where. 


1.  Panel  deflection  at  x/a  and  y/'o 
W  =  (Col  1  value).  qb^Z-Jo.,  D  ^ 


jn  o  o  O  o 


c 

c  ...  /'<-■?  •••  Af'iALYblS  (IF  SlHPLY-bUPPORTLD  C'KUIUTKUP  1(-  iANOWlLM.  _ 

C  fLaT  PAHFL.S  IIMI'bR  PKFSSIIKF  flP  PKt-S<;ilKF  PI.IK.  IM-PLuNb  LllAO;, 

C 

C  IHF  IMPin  DATA  Ai'lD  SILI(.Tbi>  (.DI-I.ST  An  f  S  API  Ab  I-DLLIIWS . 

(,  AU  =  ^I'VAPIA^'I  I  =  t  X-  (  I +1  Y/l  X  +  /  -DIIYX  1/ (  1-(1IIXY  n^ilYY  I 

(  Al,'  =  iriVAKIAMI  /  -  (,XY  -  MIIYv  •  I  X/ (  I-I^DXY'-MIYX  ) 

C  AKX  =  HIK.KI.IDI,  (.ni-bl  K  II  .''1  II  111!  X-DlPbCl  H.'N  bUU'.D  bl'.i.P  I'PJOKAii  hL'J 

(.  AkY  =  l(li(KL|Di.  I  III  bb  I(  I  bli  I  IM  iHi-  V-DlPbf  rilll'l  bllDMD  fl'i),.  i(,i' Ar,  At/> 

(,  alpha  =  AI  b  II  .sm;  I  (bX/l  Y  ) 

C  I'bTA  r  I!  r  alpha -I'UYX  +  <?"(,AK'-’A 

C  Dll  -  bX- Tb  -(  l-l■■DXY  -^'HYX)  ) 

(,  D?^  =  D1  1  ■  (I  Y/bX) 

C  ..  DAA  =  T>-  -(  bXY/^ 

(.  DUX  -  ("i.rx 

(.  DHY  =  I.XDrY 

L  I.AM.A  i:  (.  =  (.<Y=  (  I -.-HIXY  -'’IIYX  )/MiKl  (bXi-l  Y  ) 

(.  lA  ■  TlllAI  .iliHDbP  lib  A/l<  (Pa‘ii-1  ASPl  I  T  I'AIIIiS) 

i;  (I  =  mini  IIH-'IDH  III  PHI,  l/HbP(  PHI  =  I'lY/i-'X 

(,  IP  =  'DIAI  HH.il|)l  •<  (lb  V,  WHbPI  P  =  P  1 /WSDK  I  (  D1  I  -llXX  )  /  (  b-'  V’'H"X  ) 

{,  b'l  =  Tl'IMIAllllbl  lliTrl.iP  Fhk  SI-RlbS  'iHriH  Al  I(  il  i 

(,  PSl  =  S'iPT(FX=/RY)/(  I-HDXY'-MiiYX) 

K  =  DOX/DHY  =  (,CX/(.CY 
HX  =  I'X/HXCK 
PY  WY/MYCP 

XI  =  SDPTIDl  '.pD>!2)  .  - - - 

D  IMbliSIllli  ADPC  ?0)  ,  XUA(XD)  ,  Y(M(?()  )  ,  PHI  (  20  I  ,  V  (  -h)  )  ,  OX  (  )  , I’V  I  XH  )  , 

1  ADX  (  70  I  ,  ADY(  ?(■'  I  ,AMXY(20)  ,"(20)  ,AOX  (2(|)  ,aOY  (  2D  ).  SY  I  2D  )  t  (.  Y  (  20  ) 

7, APS ( (20) ,TIT( 181 ,AAKX{20l ,AAKY(?0  ) 

1  P,bAl>(  A,  n  (1  I  Tt  I  ),  I  =1 ,1  Hi 

7  Fimo/iKlPAA)  - 

"pbAiK  '5,2  )  TMAOH,  ALF,  TMI- ,  TNV,  A  I  1  ,  A|  2 ,  P  X  ,R  Y  ,  T  NPS  1  , 1  NXI  i  A  ,  IMYDP 
7  FnRMAT(6bl2.0| 

A  =  ALF 

111  SO.01AAA7A 
I  ArTOADn 

MA=TMXDA 

MP=TbiY()H 

3  RFAn(5,2)  (ADR!  I  ),  AAKX(  DtAAKYdl,  1  =  1,  lA) 

6  IPrTliv 

RFAIM  ‘i,2)  (  V(  I  ) ,  1  =  1,  IP) 

L?=TMPbI 

MFAH(R,2)  (APSK  I  ),  1=1, L2) 

PFAIMH  ,2  )  (  XHA(  I  )  ,  I  =1  ,MA  I  ...  ...  _ 

. Ri-Aii(  8,7  )  (  Y'lP  (  I)  ,  I  =  1  ,Kh  I 

Dll  100  L3=1,L2 
DD  10(1  .J=l,IA 
A«l  =A(IC.(.I  ) 

AK  X  =  AAk,X  (  J  ) 

_ ;_AKY  =  AAKY(  J  )  _  _  _ _  _ _ _ _ 

AH7  =  AHr-'*2 
AR3  =  ARl=’il  , 

AR4  =  ARl-ii  A 
PSI=AP.S1  (L3) 

(.AbP’A=(AI  l/(7.=-PSI  I  +  A12/PSI  -.‘'iHA  -.8/A) 

HFTA=( 1.8- AI l/Pbl  +  2.=AI7/Pbl  -1.8-A  -l.S/A)  -  - 

H=HFTA 
firOA.MI- A 
W I TF (h,«) 

R  FORMAT! • I ' ) 
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wniTKft,?)  (T n  (  n,  1  =  1,  iH) 

10  WRITK6,ll  ) 

tl  PlIKOAII'O  nRRL(0.  r  IflM,MIIMt-NT  ANO  SHFAft  (.OFfUf.  IRNTS  FItk  SIMPLY  SU 
KM’dinill  MU  IiiiilPIII' I(  SAi'lllW  K.M  PLAflS'l 
1  *>  HP  I  !!■  (6, 1A1P,(.,|1,A,  PSI 

I A  f  IIPMAI  (  *(11  ‘,y  •  U  X/i,<  i(.Y='  I  S./'.bX  •OAMMA='  't)HA=  '  FM.  },SX  '  ALP 

1  HA  =  'l-A.  t,  '>y  HI',  I  =  ipi  I  .!)) 

17  '-'P  I'll- ( A  ,  I  (I )  Ai|P(  J  )  ,i{X,KY,AKX,AKY 

IH  PtiPOAl  ( 'iit^x'/v/rl  ='P5.P,bX'NX/NXf.P  =  •  H  A  .  <  ,  SX  •HY/M  YC  k  ='PA.1,SXiKX 

l  =  'P‘'.?,>,yiKV  sipA.p) 

?S  0(1  (00  1  =  1  ,  |P 
?A  HP  n  I  (  6 , 7  7  I  V  (  !  ) 

71  pOP'AT  (  '  O  ipx  '  .  .  ,I)A  r  A  KIP  V  =*PH./.| 
i)M  1  00  .1  1  -  1 ,0A 
(  1  =■<.  lAlA  XllAf.ll  ) 

?S0  i/p  ni  (A  ,/’Sl  )  XIIA(  J  1  ) 

?S)  (  nP'‘'AT  (?‘.M  '■.OhPP  K.IPMTS  Kll'  X/A=F6.3  ) 

SOS  I'P  III  (o  ,  SOA  ) 

SOA  KipOATI'O  li»X  1 l'X/<i-/i,v./p  MY/0~iiv-V  h 

lXY/“-i-‘>7  ox/owl;  OY/(j;l<') 
on  I'll, 
t/(,l?)  =  .0 
AMX(,JP  )  =  .() 

AMYI J?)=.0 
AMXY(  .)2)  =  .0 
AOX(  J?)  =  .() 

AOY(.I?)s.O 
lOO  COM  I  1001- 

30  on  "o  (1=1, MI, p 
on  <>o  .'-i=I,M,P 

OMrM 

(iO->' 

iSP  MO  (0,1  1,1  I 
L  '  III'.  (Ill  III 
(III  poll  K  r  I  ,  op 

SY(I<)-'.|M  (YOPIKI  IM'ilA  OF'l 
I  Y  (  K  I  •  ro'.  (  YOH(  K  I  3  .  I  A  I  A  OM  I 
'0(1  I  o,  I  (illll 
|I^■P=I|■  / 

-j-i 

O■•’A=ll0p•.  -P 
li.'iprOr',.  =  !P 
(Pl3=ni'H:IIMp 
OOAllpiP  -.-p 

All=  (A  llO'l/AI-.P  +0'''-MPYP.  (,/AOl  )-(-l.  ) 

A1P=  1.  +  r',i',p-.-A=.  V  (  I ) /..I'P  +oi;p  •r."\/(  I ) 

A1  3  =  (-l.  !  Klin.  OiY  -VI  I  ):.-K  Kl'-(,  l/AHl  ) 

API  =  (0113/ A  +OMv(l,,-..  ,;r  ) 

AP?=  A1  3/p 

A  73=  1,+P  >-(ON7/A  +(,-.niMp/A(<P  )!•  Vdl 

A  3  1  =(  -  I  .  )  "(  A,  ll'.A/Ai'A  +P  .  ~  1'.'  Ilf  p-  (IMP/ AO/  +  nOA/AI 
A3?=  -()(,  /Alil-  (ll.iP'iA/Ai'P  <■((■"/•  i>  )■•<''(  n 
A3  3=  I'll'C  P'‘M(  I  )•' ((I.JP/A  +  1)1  P  vn//vi(P  ) 

a:  M  =  (  AlP  IAP3-APP-A13  )  /  (  Ai3v(  A  11  5A2  P-AlPv  A2  1  |  + A32“  (  A  1  3- AP  1 -i.  1 1  3  ) 

1  +  (-A31-AK.X  >l'X-  n(3//AH2-AKY=P  Y=l'M2 1  «(  A  1  Ps- A23-A/2-I A  1 .3 1  I ■■  01/lo(  s-M() 
Of!=AP2-Al3-A12=XA2  3 
ALFP=(  A11:APP-A1P*A21  )/')l) 

ALFI=(A1  1  =A?3-A13*A21  )/l)L!  _ 

S=V( I ) 

HMX=-  P.HAPA  oAMM=(  (ni3aS*ALF  1/ At3  l-iiM2/ AbP  I -» A+ (  h-2  . -G  1  “  (  K  =S» 

10M«AI.F2-nN'P  I  I 

PMY=-  S.RA9A  *AMN=(  1./A»I0M*ALF2*K*S  -()N2  )  +  ( 13-2  .  *G  I  *  (  AL  F  1  •>0M 


1  wS/AHl-lll'»w?/AH^  )  I 

KMXY=(;*  *AMN»(ON*ALf  l*S  +ALF2*R*UM*  S  -2  .  •>UM«nM/A«  1  ) _ 

HfIXr-  Sl.OOAS  *AMNv(IHvillN*  (  1  .  /  ((INeAB  1  )  a  I  ( B-2  .*0  )  *  (  UN»AL  H  2«K  "s 

l-HM?  l+Aa(ALFlanM»S/AHl  -OM2/AH2))  +r,*UNa  (  ALF  1  *b/nM  +ALF?aKvS/ I  I1N<- 
?AR1 )  -2./AH1 ) ) 

B('Y=-  il.OOAS  vAMNa  ( Ga(lM  a  ( ALF  I  aS/OM-f  ALF2  *k*  b/ ( UN«  Ali  I  ) -2  .  /  Ah  1  I 

la( l./AHl )+( ((lNaALF2*R»S-nN2)/A*(B-2.»G)»(ALFl»UM»b/ABl 
2-i'M2/Ah?  )  ) 
lit'  J?  =  1,11(5 

A.‘’'X(J2)=A(.X(J2)+t*'-'XabX»bY(J2)  .  .... 

A(iY(  J2)=A(-1Y(,I2)+I1«Y«SX=SY(.I2) 

Ai'XY  (  J?  )=Ar’XY(  ,12  )  +h(lXY>  (.X  Y  ( .12) 

A(1X{,I2)=A<IX(,I?)+|1''X-CX»SY{  J2) 

A<1Y  ( ,12  )  =AOY(  ,12  )+h''YaSXaCY  ( ,12  ) 

.•'(,I2)=1!(,I2)+A.Y('»  SXabY(  J?  ) 

100 

(If,  l=(l 

IF  ((ifil-(r(i)G,)  ,AO(),<IO 

(. 

C. ...  ,1(1111’  Til  l’|(K  IIP  ACM  tL-K.', 

1, 

400  |F(  '-0)4(11  ,4(T3, 401 

401  lli,zf< 

||M  =  H 

(,('  TO 

4nx  r.diniiiiih 
SOO  OP  I  TF  (  A,OI)|  I  0,1', 
sni  FnP'XAT(4M0  0  =  ,  I2,XX,2HO=,  1  2) 
bio  IH  b2n  ,12  =  1, '■•h 
oil  Ilk  I  TF(6,bl2)  Yoni  J2  ),W(,)2  )  ,  AhX  (  J2  ) ,  AMY1J2  )  ,  AMXY  (  J2  )  ,AUX  (  J2)  ,AOY(  J2) 
■112  F(iM.iAT(bH0Y/(l=F6.X,  F9 .  A  ,  F 1 3 , 3 ,  F  14  ,  F  14 . 3 ,  H  J  .  3 ,  H  1  .  1 ) 

b20  CnOTINUF 
90  CnOTIMUF 

1.0(I_C.()MTIflUlh.  _  _ _  _ _ 

CO  TO  1 

FA'O 


wiiere 


“H  •  *f  =^/(2  (1  - 


"22 '  “uV\ 


<l  “  pressure 


2.  Panel  moment  at  x/a  and  y/h 

^x  ~  ^  value)  .  q  1> 

3.  Panel  moment  at  x/a  and  y/b 

f'Jy  =  (Col  3  value) .  q 

4.  Panel  moment  at  x/a  and  y/b 

"  (Col  4  value) .  q  b2 

5.  Panel  sliear  at  x/a  and  y/b 

*  (Col  5  value) .  q  b 

6.  Panel  shear  at  x/a  and  y/b 

Q  =  (Col  6  value),  q  b 


EXAMPLE  PRnRT.KVt 


lioneyOTi*''^el  wither. ^0  5™  foA«%f  T  * 

leaded  under  prassu«  enly'and  tlin  !  f.o1oVN  /T"  “ 

^V-V  -  0.5.  y  X  x'^'x.cr 


y  y,cr 


Many  of  the  parameters  on  tJie  data  sheets  shovvn  are  taken  fmm  \r  c 
example  problem.  Tte  results  are  tabulated  as  Mta.r 


I 


.'MOJJ'*, 


NUMBER  IDENTIFICATION  DESCRIPTION  DO  NOT  KEY  PUNCH 


(O/^-AS/O)  S  LAMINAlh  AT  kfltlM  TFMPEKATUkF  ,  IIMOFR  PKFSSURt  ONLY 


T 
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cni-Pt-ir, ifcMTS  f-op  x/A=  ().'>r)ri 


COEFFICIENTS  FW  X/A»  0.500 


. . . . . . . 


OFHLFCT  |ll^  .HOMtWT  AIM)  <;hFAH  f.tifc' FF  I  C  I  FNT  S  HIM  SIMPLY  SUPPDRTFl)  flk  THOT  POP  i  C  SA.vllJWlCH  PLATbS 


c^l:HPI^.l^NTS  FMH  x/A=  fi.'jnn 


. . . . I*  '  '  .'W"  . . . . . 


CriEFUCItNTS  F(IH  X/Aa  O.bOO 


T 
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(0/4—4S/0)S  tAMINATE  AT  RtinH  TEMPeRATUHFt  UMOER  PRESSURE  PLUS  NX  AND  NY 

DEFLfcCTJUN.MUMENT  AND  SHEAR  r.Qfi FF IC IfcNTS  FUR  SIPPLV  SUPPURTEIJ  ORTHUTROPIC  'SANDWICH  PLATES 


T 


OCrf_  PR3GRAM  AC- 11 

H/G  SANDWICH  PANEL  STABILITY  FOR  IN-PLANE  SHEAR  LOADING 


ADVANCED  FILANe'ITARY  COMPOSITE  MECHANICAL  PIOPERTIES  COMPUTER 

PROGRAM  D1-:CK  NO.  AC- 11 


ABSTRACT 

This  program  calculates  the  general  buckling  load  for  a  flat  honeycomb 
(H/C)  panel  with  orthotropic  face  sheets  and  core  under  in-plane  shear.  Ihe 
program  listing  includes  a  modified  eigenvalue  solution  subroutine  whicli 
involves  convergence  of  the  geometric  mean  of  succeeding  eigenvalues.  The 
face-sheet  thicknesses  are  assumed  to  be  equal. 


PROGRAM  LNPUI'  DATA 

The  input  data  can  be  of  two  types,  material  property  input  or  material 
parameter  input.  The  input  data  format  is  all  6F12.0  except  for  the  title 
card.  Sample  data  sheets  are  included  following  the  text  in  this  section. 
Tlie  data  description  is  as  follows: 

Card  A: 

TIT  =  any  desired  one-card  title.  The  title  is  in  columns  1-72. 

If  no  title  is  desired,  a  blank  card  must  be  used. 

Card  B: 

CNTL  =  Input  data  fcrm  indicator  which  is  0.0  for  material  property 
input  (cards  Cq  and  Do)  and  1.0  for  material  parameter  input 
(cards  Cj  and  Dj). 

Q  =  Range  of  m  and  n  (maximum  of  13.0). 

Card  Cq’.  Used  if  CNTL  =  0.0. 

EX  =  Young's  modulus  in  the  x-direction. 

BY  =  Young's  modulus  in  the  y-direction. 

GXY  =  Shear  modulus. 

UMXY  =  Poisson's  ratio  v  . 

xy 

GCX  =  Core  shear  modulus  in  the  x-direction 
GCY  =  Core  shear  modulus  in  the  y-direction  (G'  ). 


Card  D^:  Used  if  CJTL  -  0.0. 

A  =  Length  of  panel  in  the  x-direction,  a. 

B  =  Width  of  ptinel  in  tlie  /-direction,  b. 

C  =  Core  depth,  c. 

TF  =  One  face-slieet  thickness,  t^. 

XNIT  =  Maximum  number  of  iterations  to  be  allowed  (usually  set  at 

100.0). 

Card  C,:  Used  if  CN'TL  =  1.0. 

.•\LPiiA  =  The  alpha  value  from  program  AC- 2  ^ 

PSI  =  Modulus  ratio  4>  found  from  AC-2  where  ip  =  \ /b  b”  ! (l  -v  u  ) 

V  X  y  /  '  xy  yxi 

V  =  Core  shear  parameters  where  V  =  +  t^) ^(2G(  b“j. 

R  =  Core  shear  modulus  ratio  =  G'  /G'  . 

cx  c>' 

All  =  Invariant  1  from  AC-2.  I ,  =  b  /l  +  b  /b  +  2v  v  \. 

1  x\  y  X  yx//\  x>'  yy.) 

AI2  =  Invariant  2  from  AC-2.  I-,  =  G  -  \l  v  /( 1  -  u  u  ). 

xy  X  yx/l  x>'  yx  / 

Card  U^:  Used  if  CNTL  =  1.0. 

A0B  =  Panel  iispect  ratio  a/b. 

C  =  Core  depth,  c. 

TF  =  Tliickness  of  one  face  sheet,  t^. 

XNIT  =  Maximum  number  of  iterations  to  be  allowed  (usually  sot  at  100) 


182 


Hie  input  data  floiv  chart  for  AC- 11  is: 


CASE  CASE 


Additional  cases  can  be  stacked  one  liehind  the  other. 


OUTPUT 


The  program  output  is  of  two  t>'pes,  depending  on  uiietlier  material  property 
or  material  parameter  input  is  desired.  For  material  property  input,  tlie 
output  includes  the  input  data;  tlie  material  parameters  a,  tp,  and  V  resulting 
from  these  data;.  tJie  number  of  iterations  (NIT)  carried  out  for  ks^gven 
kg  and  tlie  final  buckling  load  per  inch,  For  material  parameter 

input,  tlie  input  data  are  printed,  as  well  as  the  two  values  of  NTT  mentioned 
previously,  and  the  final  buckling  load  as  a  function  of  l/b^,  since  for  this 
input  no  b  was  introduced. 


c 

C  PROKR/H  ACll  ...  H/C.  P6NPL  bTAHILtlV  FOR  IN-PLAIlfc  ShFAR  LOADING 

C 

C  THF  I<MPin  DATA  AND  iK^CTU)  CtlNSTAMTS  AR6  AS  FOLLOWS..... 

c  A  =  pamil  widih  along  X-DIRFCHON 

C  AH  =  IwvAKIANI  1  =  FXsn+6V/FX*2«MOYX)/(  l-MOXY*MOYX) 

C  AI/>  =  iMVAKlAt-lT  ^  =  (,XY  -  MOYX^FX/ (  l-HOX  Y«HO  YX )  _ 

C  ALPHA  =  S('KT(FX/FY) 

C  Aod  =  PAtIrL  ASPFCI  KATID  A/B 

C  h  =  PANFL  LFHMH  ALONG  Y-f)  IKFf.  T  ION 

C  C  =  CORF  depth 

C  CNTL  =  f.ONTRIH.  card  FOK  OF S I GNA T  ING  TYPfc  OF  INPUT  T(i  b6  USED 

C  1=0  FOR  MATL.  PROPtRTIfS,  OR  =1  FOR  .iATL.  PAKAHeiERii- 

C  FIi;mio  =  Fir.FHvALOfc  ITERATION  SOLUTION  THAT  TAKES  THE  SORT 

c  OF  T‘iE  Ahs  valuf  of  soccEfniNG  figenvalue; 

C  GCX  =  CORF  SHEAR  MODULUS  IN  ToF  X-DIRECTION 

C  GCY  =  ccRF  Shear  moooi.us  In  The  y-oirection 

C  PSI  =  SORT(EX=EYI/tl-MOXYoMOYX» 

C  O  =  NdhbFR  of  TERMS  IN  ThF  <;FRIFS  IMAXIMUH  of  13) 

c  R  =  Gt.x/i;rY 

C  TF  =  imr.KNESS  OF  ONI  FACF  SHFFT 

C  TITII:  =  ANY  OFblRFO  ONf-LiNF  TITLE  IN  COLS. 1-72 

C  V  =  P|**2*PSI«TFa(C+TF)/(?*GCX»B**2» 

C  XOn  =  MIIMHFU  OF  ITERATIONS  OFSIREO  IN  THE  SOLUTION 

r. 

OIMF'ISION  H(  MS).  TITIF(IB),  PBI2) 

10  WRl  U  (h,‘’(l) 

50  FORMAT ( iHl ) 

RFAOIS.IOO)  title 
100  FORMAT (IMAA) 


hrufia.iiio)  title 

C . RFAD  material  PROPEKTIFS 

RFAOIS, 11(11  CNTL.O 
llO  FHRKAT{6E12.0) 

IFICNTL)  101,10V, 101 
C.....RFAn  PARAMETRIC  PPOPETIES 


101  REAUIS.llO)  ALPHA, PSI, V,K,A11,A12,A0B,C.TF,XNII  ..  _ 

Cr,l=AI  1/PSI-AlPHA-l  .0/ALPHA 
BFTA  =  1.5*Cr.l*2.«AI2/PSl 
GAMA=0.5»CC1+AI2/PSI 
0MXY=ALPHAa(BFTA-2.0*GAMA) 

■jMYX=UMXV/ALPHA3a2 
AM0A  =  1  .0-UMX  YailKYX 
EX=ALPHASPSI cAHOA 

EY=FX/ALPHAS=7  _ 

f.X  Y  =  A  I  2+0'1YXCFX/AMnA 
SX  =  V*C/  IC  +  IF  !/ AlmCV 
WPlTF(A,H‘i|  alpha, PSI, V,M 

115  FDMmatI 1H0AX,2RHTHF  material  PARAMETERS  ARE , // lOX ,6HALPHA=F6. 3, 5X, 
I AMP5I=F12.5,5X,2HV=F7.5,5X,HHGCX/GCY=F6.3) 

•wP  I  TF(E,,  UN)  AOjj.CtTF  _  _  _ _ 

115  EORMATI 1hO,AX,2hhTHE  GEOMETERY  PROPER! I ES  ARE , //12xV4HA/B=F7 . 3,6x1^ 
l7Hr.  =  F7.5, 10X,3HTF=F7.4t 
GO  TO  131 

C . RFAO  MatfhIAL  PKOPFTIES 

lOO  RFAn(5,110l  FX,EY,GXY,UMXY,GCX,6CY 

C . RFAO  PANEL  GFOHETRY 

II?  REA()(5,110)  A,B,C,TF,XNIT 
Ar'H  =  A/B 

OMYX=EY/FX*UMXY 
AIPhA=(FX/FY )*e.5 
AMOA=) .0-OMYX*OMXY 
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PSl=(FX*FY»#«.S/AMnA  _ 

(;ama=(;xy/psi 

HF  T  A=ALPHA*I>MYX*?  .n*r.AMA 
P  =  r.*TF 

1=PS1«TF 

l)= 

SX=C.ini--.Hn‘>#;/ (PcAltS^/r.CX 

R=(;f.x/(;r.Y 

V  =  SX9psAllil-;*7/r. 

C . WPIIF  ifJUilT  riAIA 

118  W8  I  ff  (h,  IPO)  FX,(;f.X,FY,(;CY,GXY,UMXY,l)MYX 

1?0  FOKMATI lHO,4X,PHHTHfc  MATgHIAL  P80PFRT16S  ARE  //lOX, 15HFACh  SMtFT  fc 
IX  =Pn.S,10X,  ‘.oCHfc  GX  =C13.5/PlX,4HfcY  =ei3.;),)L5X,4HGY  =cl3.5/20X 
?SH(  XY  =Fl3.‘i/lQX,.->HHllXY  =FR  .4/ 1  9X  ,4HH«JYX  =FK.4) 

C 

HPITF(6,130)  A.btC,  TF.AOt) 

130  FORMAT! 1H0,4X,P2HTh6  PAN.  L  GFOMfeTRY  IS  // 1  OX , 16MPANgL  LENGTH  A  =Fb 
1 .3,SX, 15HPAMFL  WIDTH  H  =F«. 3/12X, lAHCflRE  DEPTH  C  =FH.3,  SX,25HFACE 
Z  SHEET  THICKNESS  IF  =rB, 3/21Xt5HA/B  =P8.3) 

WRITF(6,125I  alpha, PSI.VfR 

125  FORMAT! lH04X,2flHTHF  MATERIAL  PARAMFTFBS  ARE , //lOX ,6HALPHAsF6. 3, 5X, 
14HPSI=£12.5,5X,2HV=F7.4,5X,8HGCX/GCY,  F6.3) 

C . CALCDLATF  CONSTANTS 

131  SY=RCSX 
P=C*TF 
TrPSlsTr 
n=T9P=t=?/^,0 
A082  =  A(ili--«2 
A(rR4iA(lb29«2 
I0=ii+,00l 

C . XLMO  constants 

XLMN)  =f,AMA9Ann2 
XL'"i29  rt  /ALOHA^AObZ 
XI  nm=  R  »r,AMA 

C . XKMM  CONSTANTS 

XK«-''l  =2.0*HFTA*A08? 

XKVM2  =  A(li<4/ALPHA 

c . fmn  constant 

XFMNl=l .0-OETA«»2 
c...,.vMo  constants 

XVMfil=Atlb2*SX 

XVMN2=SX9SY 

C . DF  roOSIANT 

0FrTF993/17.0*PSI 

oFoo=nF/n 

C . Calciilatf  vAhiahLf  constants  after  Checking  !m*n( 

IPrl 
134  1=0 

NIlrXNl T 

140  fin  200  M=1,I0 
AlsM 

MA= ^.?42^7oAl 
A2  =  A  >  092 

A4  =  A?<-»2 

!F!0-?t!H/2n  145,150,145 
145  NS=1P 

_  GO  TO  155  .  . 

150  NS=3-IP 
155  on  200  :i=NS,I0,2 
l  =  I-‘l 

B1*N 
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_ _  _  _  _ 

B4bB2**2'~~ 

C . CALCtlLATE  VARIABLE  CONSTANTS 

XLMN  >ALPHA*A2'»XLNN1*B2'»XLNN2«B2«XLNN3*A2 
XKMN  sALPHA«A4'»XKMN1«A2*R2«XKNN2*94 
FHN  =  XFMNt«GANA«XKHN 

VHN  =  (  XKHN-»  ( XVNN 1  *B2'»SV*A2I*FMN }/  a.0«XLNN*SX«XVNN2»'MN ) 

vmn-!  =nFnn*xKMN 

_ AA;;<2.0»VF.NI>VHNI/A0H  _ _ _ 

c . check  to  see  if  (m«ki  and  in^si  are'ooo  or  even 

J=0 

■  on  190  lR=i.io 
RI^IM 

_ RI2=RIe«2  _  ..  .  _ 

IF(IR-2»«iR/2M  160,170,l<i0 
160  1SS=IP 

f.n  Tn  mo 
170  ISS=3-IP 
180  on  19(<  IS  =  ISS.IUt2 

_  SI=!S  „  .  _ _ 

■hC=mh»S1 

SI2=S1=”2 

J=J*1 

lFCK*I«-?»nH*TR|/2l  }  182*187.182 
182  1F<N*IS-?»UN+IS)/2I1  186*187.185 

_ 185  M11,J»=-mC  /<A2-R121/IB2-S12}/AA _ _ 

~  Gn  TO  190 
187  HI1,JI=0.0 
190  CnuTl-JOF 
200  cnmiNtiE 

102=1 l0»*2>2-IPl/2 

_ ^OFLTAs.OOI  _ _ _ _  _ 

call  EIf.HlM102.0flTA*H,e8IIPl*NlTl 
HRnFl6,2171  NIT 
217  format aH0*5X*<.HNlT *16/1 
IF.l  ip-21220, 240*240. 

220  IP=2 

_ on  TO  135  . .  .  _ _  .  .  _ 

C*..*.PP1>'T  OUT  EVEN  ANO  OOP  KS  HHICH  ARE  BBIll  AND  BBI2I 
240  BPll 1=88111/6082 
B8121=8B12I/A0B2 
IF|RHlll-8B12n  245*250*250 
245  XKS=RP.Ill 

_ GO  TO  25.5 _ _  _ 

250  XKS*8B;2I 

255  IFICNTLI  256*260*256 

256  XKNXCR>XKS«9,8696«0 
WRITE 16 *258 I  XKNXCR 

258  FnRNATl»0*5X»NSCR=KS  •  Pl**2  •  Ml  •  TF  •  IC*TF)**2/|2  •  8**2I  =• 
IE13.5*/B**2*  1 
GO  TO  280 

_ 280  XNCRS=XKS*9.R6R6»0/»**2  . 

WRITE16.2651  10*88111*88121 

265  FnRMATIIH0*38H**«*THE  VALUES  FOR  KS  WHEN  THERE  ARE  I3.22H  CNN  COE 
IFFICIFNTS  ARF,/20X, lOH  KS  EVEN  *E13.5*10X*  9H  KS  UPO  *E13.51 
HRITF16.270)  XNCRS 

270  FnRMA7l»0*5X*NSCR=KS  •  Pl**2  •  PSI  •  TF  •  1C4TF|**2/12  •  H9*2I  •• 
II  13.51 
?8U  «0  TO  10 

ENII 
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SUAH nil  TINE  EIGMlN(NNfOELTA,BlNV,AUAMIN,NlTI 
DIMENSION  C(  R5),fll(  85)<BINV(  85tA9> 


I  1  =  1 

1  nn  ?  i=i,NN 

C  (11  =  1.0 

2  CONTINUE 
N=1 

_  3  on  S  I  =  1 , NN 
■■■■81(11=0.0 
A  nn  s  j=i,NN 

81(11=81(1 )+81NV( I .J1*C( J) 
5  CONTINUE 

C1=ABS(B1(1)) 

_  ALAHnA=Cl  . . 

s  on  A  1  =  1  ,nn  . . 

c(n=Hi(i)/ci 

9  CONTINUE 

IF(N-l)  llilO.ll 

10  N=N+1 

_ BLAMnA=A|.AMnA 

CA=(HLAMnA>»'LAH0A)**.5  . 
f.O  TO  3 

11  C2=(ALAMnA«RLAH0A )**,5 
C3=A8S(C?-C4) 

IF(C3-0ELTA)  14,14,12 

J,2_BI  AMnA=ALAMnA_..  _ _  . 

C4=r.2 

11=11+1 

IFdl-NlTI  13,14,14 

13  GO  TO  3 

14  NIT=II 

_ ALAHlNal.O/,  ' _ 

■rf'turn"’" 

END 


-a 
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NUMBER  IDENTIFICATION  DESCRIPTION  DO  NOT  KEY  PUNCH 


T 


EXAMPLE  PROBLEM 

'  Consider  a  honeycomb  panel  to  be  made  of  [0/*45/0]g  face  sheets  with  a 

core  depth  of  1/2  inch,  G'  =  126,000  psi,  G’  /G*  =2.0,  and  an  overall 

cx  cx  cy 

panel  size  of  15.0  x  15.0  inches.  To  illustrate  the  use  of  the  program,  both  ^ 

I  types  of  input  will  be  used.  The  material  parameters  and  material  properties 

are  taken  from  the  AC- 2  and  AC- 5  sample  problems. 

The  data  sheets  are  shown  for  both  types  of  input  near  the  end  of  this  f 

section.  Note  that  for  the  material  parameter  input,  3 

N  =  5.1021  X  10^  /  b^  ; 

s,cr 

and  if  b  =  15,  as  in  this  example,  4 

N  =  5.1021  X  10^  /  (15)^  =  22676  Ib/in.  1 

s,cr  2 

as  found  from  the  output.  | 

§ 

I 
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NUMBER  IDENTIFICATION  DESCRIPTION  DO  NOT  KEY  PUNCH 


(0/+-45/0)S  LAMINATF  at  room  TFMPgRATURe 


THE  HATERIA!.  PROPERTIES  ARE 


FACE  SHFFf  FX 
FY 
r,XY 
MUXY 
MUYX 


O.lhhTOE  08 
O.'SnSHOF  07 
().4^0n0E  07 
0.701'; 
0.21PR 


THE  PANEL  GEDHFTRY  IS 


CURE  r.X  =  0.12600E  06 

GY  =  0.63000E  OS 


PANFL  LENGTH  A 
CORE  DEPTH  C 
A/O 


15.000  PANEL  WIDTH  h  =  15.000 

0.500  FACE  SHEFT  THICKNESS  TF  =  0.042 

1.000 


THE  MATERIAL  PARAMETERS  AREt 


ALPHAS  l.fllS  PS1=  0.10794F  08  V=  0.0423  &CX/GCY=  2.000 

NIT=  7 


NIT=  11 


#*«*THE  VALUFS  for  KS  when  there  are  13  CMN  COEFFICIENTS  ARE. 

_ _ KS  EVEN  =  0.7R484E  01 _ _  KS  ODD  =  p.fl4673E,  01. 

NSCfl=KS  *  Pl**2  *  PSI  *  TF  *  (C+TFI**2/(2  *  B»*2)  =  0.22670E  05 


0/+-45/0)S  laminate  at  room  TEMPERATURE 
THE  MATERIAL  PARAMETERS  ARE. 

ALPHAS  1.‘R15’“'  PSI=‘0. 107<J4E  08”'  V=  0.0423  GCX/GCY=  2.000‘ 

The  GEOMETERY  PROPERTIES  ARE 

A/B=  1.000  C=  0.5000  TF=  0.0416 

NIT= . ~7  ' 

NIT=  11 

NSCRiKSr*  PI**2''*“PTr  *  TF'O'C'VtF  I442/I2  •  B»*2)  '«  ^O.SloYlE  07/84*f” 


rn'/fPIlTF.R  PROGRAM  AH- 17. 


HAT-  AND  Z-STIFFENED  PANEL  STABILI'rY 
AND  STIFFNESS 


AC-12 


ADVANCED  FILAMENmY  COMPOSITE  MI^QIANICAL  PROPERTIES 
COMPUTER  PROGRAM  DECK  NO.  AC- 12 


ABSmCT 


Conputer  program  AC- 12  is  designed  to  provide  tlie  optimum  minimum  length 
and  tlie  resulting  face  slieet  and  stringer  stresses  for  a  hat-  or  Z-stiffened 
panel.  These  quantities  are  given  as  a  function  of  load  level  and  stiffener 
spacing  as  well  as  otlier  geometrical  and  material  properties.  Ihc  loading  is 
uniaxial  in  tlie  x-direction. 

Tile  optimum  minimum  lengtlis  of  tlie  panels  are  found  by  equating  tlie  local 
and  Euler  buckling  equations  and  solving  for  tlie  lengtli.  A  program  listing  is 
given  following  tlie  text  in  tliis  section. 


INPUT  DATA 


The  input  data,  as  shown  on  the  sanple  data  sheets  in  this  section,  arc 
given  more  conpletely  in  the  following  listing: 

Card  A: 

BMAT  =  Any  desired  one-card  alphanumeric  title  limited  to  columns  1-72. 

Card  B:  Face-sheet  properties 

All  =  E  t/Zl- V  V  ]. 

X  '\  xy  yx/ 


A16  =  Defined  by  equation  9  of  this  volume  and  is  0.0 
symmetric  laminates. 

A22  =  E  t/(l-v  p  ]  =  /E  /E^)  A... 
y  xy  yx/  \  y  x/  11 

A26  =  Defined  by  equation  9  of  this  volume  and  is  0.0  for 
symmetric  laminates. 

A66  =  G  t 
xy 

Card  C: 

Dll  =  Flexural  rigidity  of  the  sheet  in  the  x-direction  (load 
direction) . 

AEXF  =  E^  =  Face-sheet  modulus  in  the  x-direction. 

X 
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T  =  Face-sheet  thickness. 


TP  =  Indicator  for  Z  (TP  -  0.0)  and  hat  stiffeners. 

AK7  =  Buckling  coefficient  from  equation  131  of  this  volume. 

AKS  =  Stiffener  buckling  coefficient  found  from  ecjuation  132  of 
this  volume. 

Card  D: 

AKF  =  Face-sheet  buckling  coefficient  defined  by  equation  132  of 
tills  volume. 

Card  E: 

AEXS  =  Young’s  modulus  of  the  stiffener  in  the  x-direction. 

TNBS  =  Total  number  of  stiffener  spacing  to  be  considered. 

TNNX  =  Total  number  of  N  loadings  to  be  considered. 

DS0HS  =  Hat  geometric  parameter  dg/hg  (0.0  for  Z-stiffener,  i.e., 
when  TP  =  0.0),  as  defined  by  the  following  sketch: 


Card(s)  F: 

XBS(I)  =  Stiffener  spacings  (bg)-  (maximum  of  30). 

Card(s)  G: 

BNX(I)  =  Load  levels  to  be  considered,  N  (maximum  of  30). 

i 
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c-ionc'.  rinr>oooooooo<-><->r>or>oor>rir> 


AC12...Z-  AND  HAT-STIFFENER  PANEL  STABILITY  AMO  STIFFNESS 
the  INPUT  DATA  IS  AS  FOLLOWS . 

tit  :  "xyIt"’  ""  ^"''^TION  <,  ,0.0  POP  SYMMETRIC,  LAMINATES) 
T  =  LAMINATE  THICKNESS 

AEXF  =  face  SHEFT  YUUNG'S  MniiULtlS  IN  THE  X -I)  I  RE  (  T  t 
TP  -  INOICATUK  EUR  7,TP=0.0)  AND  HAT . TP= 1 ' 

ARC  ^  Cf'ff'PfCIENT  FOUND  ERUM  FOUATIUN  ITl 

ARC  °  buckling  COEFFICIFNT  found  FROM  FOUATION  132 

Acyc"  CUEFFICIENT  FOUND  FROM  FOUAT  ON  32 

AFXS  =  YOUNG'S  MODULUS  FX  OF  THE  STIFFENER 

TNBS  =  TOTAL  NUMBER  OF  STIFFENER  SPACINGS 

0'!nHC^'HIT^??cc^''  '-OAOINGS  TO  BE  CONSIDERED 
STIFFENER  GEOMETRY  RATIO 
BNXI I)  =  LOADINGS 
XBS(I)  =  STIFFENER  SPACINGS 

dimension  BMAT 1 1 fi ),BNX( 30 1 .XBS(30  I 

1  RFAD(5.2)  ( BMATIMI , M=1 , 10) 

2  FORMAT! IBAA) 

A5  WRI TE (6, A6) (BMATIMI ,M=l,lfll 
A6  FORMATIlHl tlRAA) 

RBAO  IN  STIFFFNER  PROPERTIES 
5  REAn(5,3)  AEXS,TNBS,TNNX,0S0HS 

7  lAaTNBS 

IB«TNNX 

ReAD(5.3l(XBS{I),I=l,lA) 

ReAOI5t3)IBNX(I),l=i,iB) 
on  200  J=1,IR 
SO  WR1TE(6,S1 IBNXIJi ,T 

II  SRIT?j6!tn  NX(LB/INCH)=F7.1,7H  T  =  F7.A, 

Al  FnRMAT(99H0  BS  HS  TS  FS  D)  .  , 

IRIT.  SIGS  SIGF  DS 

00  200  1  =  1, lA 
ANX=BNX(J I 
BS  =  xeS( I  I 

BETA=  (AK7/AKS)«*,5 

ALFA  =(AFXS*AKF/(AFXF*AKS1)»*.5 

/F<TP)17,16,  17 


T=F7.A) 


local  buckling 

CALCULATION  OF  GEOMETRY  OF  Z-STIFFENER  ELEMENTS  TO  INSURE 
simultaneous  LOCAL  BUCKLING 


198 


o  o  uo  o  oo 
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16  MS  =((  ANX-fAfiXS/l  AKS*(ALFA«T/RS)**2)  -AEXF*T )  *nS«*2/ ( ALFA*T  *AEXS 
1*(2.*RETA  +1.  ) )  )*«.? 

TS=(ALFA  a  MS  *  T)/PS 
FS=  P.ETA  *HS 
r,n  TO  20 

CAI.r.llLATinN  OF  OnOMFTP.Y  OF  HAT-STl  FFEMRO  ELEMENTS  TO  INSURE 
SIHMLTAUEmiS  LOCAL  RIICKLIMO 

17  IIS  -((  AMXaAFXS/(AKS*(ALFAaT/PS)**2)  -AEXF*T )*nS«a2 / ( ALFA*!  aAEXS 

l*2.a(l.5+nsnns))  )*■-“.? 

DS=  nSOIlSaHS 
TS=  ALFA^HS  «T/BS 
FS=.‘^*  IIS 
20  IF(TI>)0,n,<} 

P  AKl=AEXSaTSa{2.*FS+IIS) 

AK2-AEXS«TSa(FSa(|iS-*-.Sa{T-TS)  )  +  .  6aMS*(  MS+T  )*FS*  .  5a  (  TS+T  )  ) 
Cl=-AFXS«TSa(  FS«(MS  +  ,5a(T-TS)  )  +  .5atiS«(HS-^T)  +  FS*.5*(  TS+T)  ) 
C2=-AEXS»TS*(F5«(MS>'<(MS  +  T  -TS)  +TSa(  TS/T.-T/2.  ) +T«*2/A .  )  +FSa(Ta 
1*2/^.  +T*TS/2.  +  TSa  =  :2/3.  )  +MSa(IISaa?/3.+HS«T/2.  +T*a2/A.|) 
on  TO  10 

9  AK.1=  2.«AEXS»TSa(nS-^HS+FS) 

AK2=AEXSaTSa(nS*(TS+T)*HSa(HS+T»  +  2.«FS  >'  (MS  +T/2.  -TS/2.)) 

Cl=-  TSaAfiXS  «(  llSa(TS+T)  +MS«(HS+T)  +2.aFSa(MS  +T/2.  -TS/2.)) 

C2=  -AEXS  az.^TS  a(  ESa(HSa(HS+T-TS)  +TS«(TS/3.  -T/2. )+Ta*2/4.  ) 

1  ♦r»Sa(,TS**2/?.  ♦TaTS/2,  -t-TaaE/A, ) 4)15* (M^*a2/3.  AHSaT/E.  ♦Taa2/4,)) 
10  nELa(AK2/nS)’ /(  (A12-Ai6aA267AMi»*a2  -(AKl/BS  4A11  -AIM*2/A66)  a 
l(A22-A26a*2/A66)) 

BETA*  nELa(A22-A26aa2/A66) 

01=  -(C1/os)*pfta-c2/bs  +  nil 
ANXL2  =  0. 0606*01 

XLCP=(AMXL2/AHX)aa.q  -  .  _  .  .  - 

OETERMIMATIOM  OF  STIFFENER  AMO  SHEET  PROPERTIES 
IF(TP)1A, 12,14 

12  C3=  AEXF*T  +  AEXSaTS* ( 2.aF S+  HS  -2.*TS)/BS 

0(1  TO  15  ■'  -  .  ■  •  - 

lA  C3=  AEXpaT  +  AEX SaTSa2,a ( nS+HS-*-  FS)/BS 
15  SIGF=  AMXaAEXF/C3 
Sir,S=  AMX»AEXS/C3 

mo  WRITF(6,101)PS,HS,TS,FS,ni,XLCP.  ,  SI  PS  ,  5  ir.F,OS 

101  FnRMAT(F6,3,3F0.3,F14,2,FlA,l,2Fll.1,F16.3)  _  .  . 

200  CONTI  HUE  ■  ■ 

GO  TO  1 
EHO 
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Uir  OIJE-LINE  TITLE  IN  COLS.  1-72 


EXAMPLE  PROBLEM 


An  example  problem  is  given  for  a  liat-stiffened  and  a  Z-stiffened 
[0/t45/0]g  face  sheet.  The  input  values  are  given  on  the  data  sheets  at  the 

end  of  this  section.  For  tlus  laminate,  from  AC- 2,  we  take  tlic  face -sheet 
properties  as: 

E  =  16.67  X  10^  psi 

Ey  =  5.06  X  10^  psi 

G  =  4.30  X  10^  psi 
xy  * 

V  -  0.7015 
xy 

V  =  /e  /E  W  -  0.2129 
yx  \  y  X/  xy 

The  laminate  thickness  is  8  plies  x  0.0052  inch/ply  =  0.0416  inch. 

Tlie  buckling  coefficients  k^,  k^,  k^,  are  found  from  equations  131  and  132  of 
this  volume. 

The  results  are  given  first  for  the  Z-stiffened  panels  and  then  for  the 
hat-stiffened  panels.  The  panel  critical  lengths  may  be  read  directly  from 
these  results. 
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CaviPUTER  PROGRAM  AC-40 


N-PLY  LAMINATE  THERMOEIASiTC  ANALYSIS 


AC“40 


ADVAJCrD  FIJ.J\MEOT;'J5Y  CCMPOSITE  ^4ECHANICAL  PROPERTIES 
COMPUTER  PROGRAM  DECK  NO.  AC-40 


ABSTRACT 


This  program  calculates  the  thermal  stresses  and  coefficients  of  thermal 
expaiision  in  both  the  individual  ply  axes  and  the  laminate  axes,  given  the 
macromechanical  properties  of  each  ply.  The  program  is  now  set  up  to  handle 
boron/epoxy  laminates  at  room  temperature  (R.T.),  220°  F,  and  350°  F.  If  other 
materials  or  temperatures  are  involved,  the  form  of  input  may  be  altered 
accordingly. 


PROGRAM  INPUT  DATA 

The  program  input  is  minimal  and  consists  of  only  two  cards  per  case, 
which  are  shown  on  the  sample  input  sheets  following  the  text.  The  data  are 
described  as  follows: 

Card  A; 

N  =  Number  of  plies  per  set  in  the  laminate. 

M  =  Material  temperature  indicator.  1  =  R.T.,  2  =  220°  F,  and 
3  =  350°  F 

DT  =  Temperature  increment. 

Card  B: 

TH(I)  =  Angle  of  ply  I.  There  are  six  of  these  to  a  card  (maximuin 
of  20). 

Note  that  only  one  "set"  of  a  laminate  need  be  coded  in.  For  example, 
for  a  [0/±60/902]3s  lan»inate,  only  the  0°,  +60°,  -60°,  90°,  and  90°  plies 
need  be  input,  giving  a  total  N  of  five  plies. 


OUTPUT 

The  output  consists  of  the  DI,  B,  a,  and  B  matrices;  the  coefficients  of 
thermal  expansion  in  the  laminate  axes  ax»  “y»  and  SxyJ  the  thermal 
stresses  in  the  laminate  axes,  Oy,  and  plus  the  thermal  stresses  in 
tbe  individual  ply  natural  axes,  o‘.p,  and 

The  output  is  shown  on  the  sample  output  pages  following  the  text. 


o  o  o  o  o  o 


PROGRAM  ACAO,,,N-PLY  LAMINATE  THERSQELAST IC_AMALyiliL.J>.BQGRAll_„.  . 


INPUT  DATA  IS  AS  FOLLOWS . 

N  =  MMMRtR  OF  PL  Its 

M  =  MATL.  PROP.  TFHP.  INI) IC ATOR . . . . M  =  1  FOR  R.T.,  M=P  FUR  220F , 
C  M-1  l-lll'  Oilll 

c  m  ■  T£Mi'r:i' Aruue  iN(.R6MENr  in  otcRuts  r 

(.  TH(|1  ■  INOIVIOUAL  PLY  ANGLFS  THgTAIl) 

I 

DIMfNSlOM  It  I  201  ,tT(20l ,6LT(201 lULTI 20  I .UTt 1 20  I , aL I 2o 1 1 aT I 20 I , 

1  HI  ?f))  ,  A(2o,3i3  I  ,SL  I2U1  tST(20>  ,SLr(2())  ,  10(20) 

OIOIHSIOH  TIPI),  3t3|  tPM(20l  ,EM(20)  ,IM20t  -itS)  ,01  (20,  1 )  ,AX(  2(1)  f 

1  AXY(2n) ,I(B( ^,3  )  ,SX(20) ,SY(20) ,SXY(20) ,BI (3,3) ,AYI 20) 

OIMFMSIOM  (1(20,3,3) 

FOni VALFMCF  (0(1), 01(1)) 

5  RFAO( 5,10)  M, M.DT 
10  Fni<MAT)21  12,  F12.0) 

C  TH)I)  =  AIKILF  BFTWFEN  NATURAL  AXES  OFI-TH  PLY  AND  LAMINaTP  AXFS 
15  RFA0(5,20)  ( TH( 1 ) , 1=1 ,N) 

20  FniiPAT(6F12.0| 

C  PKOPERTIFS  OF  IMOIVIOUAL  LAMINA 
AO  on  100  I=1,N 

C  CALCUIATF  FLFMFNTS  OF  T  MATRIX 
IP(M-2)A2,A3,AA 

C  ROliM  •ll-ltPPP AlUkF  PROPERTIES  _  _  _  _  _ 

A2  PL ( I ) =2DOOOOOO.  '  . . 

FT  (  I  )=271()00n. 
r,LT(  I  )  =700000. 

ULT(I)=.21 

UTL( I )=0LT( I )*ET(I )/£L( I ) 

AL( 1 )=.0O000?3? 

ATI  I  )=.0000106'7  . . .  ‘  ‘ 

).(  I  )=.O0‘i2 
on  TO  A6 

C  220  OFGRFF  PRIlPFRTIES 
A3  ELI  I )=29900000. 

ET(  I  )  =  1A20000. _  _  _ _  _  _ 

r.l.T(  I  )=377100. 

OLT(  I)  =  .21 

UTL(  I  )=ULT(  I)*ET(n/EL(  I  ) 

AL( I )=.0n0002AlA 
ATI  I  ) =.0000I39A 
H(I)=.0052 
on  TO  A6 

C  350  DEGRFF  PPnPFRTIES_  _ 

AA  ELI  1  )  =29900000.'’ 

ET( I )=1130000. 

GLT( I  1=300000, 

ULT(  I  )  =  .?! 

OTL( 1 )=ULT( I )*ET(I )/FL( I ) 

AL(  I  )  =  .00000?R3.  _  _ _ _ _  _ 

■ . "AT(  I  )  =.0000197R  ■■  ■■  ■  ■ 

H( I )=.0052 

A6  T(I,1,1)=FL(I  )/(l.-Ui.T(I  )*IITL(I)) 

T(l,l,2)  =  EL(l)*UTL(I)/(l.-0LT(n*UTL(l)) 

T(I,1,3)=0. 

T(I,?,1)=T(I,1,2) 

T(  I  ,2,2)=FT(I  )/ (  1.'-0LT(T  )*nTL(  I  )■) 

T( I ,2,3)=0. 

T(I,3,1)=0. 

T(I,3,2)=0. 


T( It3f3)»GLT( I) 

C  CALCULATI:  ELFMENTS  Ot_^MATKl>!  _  _ _  _ _ _ _ 

EM(  i  ).C()SD(TH(  i  n 
EM( I )=SIMO(TMt I ) ) 

50  H(  1  ,1  .U  =  1!  It  Itl  )«FH(  n**4+2.*(T(  1,1,2)+2.*T{I  ,3t3)  I  *  ( EMI  I  1  «»2  I  * 

1  (  FN  (  I  )  ■>»?  ) +T  1 1  ,2 , 2  )*EN  (  11  *«>4 

B.I,l,2):^(TII,l,l)+T(l,2t2)-4.*T(I,3,3l)*lEMn  )**2  1 « (  FNI  I  )  v42  ) +T I  I 
_ 1 ,  1 1?)*(EM(  n<"l'4+EMl  I  )**4) 

B(Itli31  =  IT(I,l,n-T(l,l,21-2.*T{I,3t31)*ENlI)*bMin*<‘3  +  (T(l,lt2)- 
1T(  I  ,2,2  )+2.4T(  1,3,3)  )#FM(  I  )*FMI  1  )*’»3 
B(l,2,2)=Tn,l,l)’>FN(I  I  ■>■>4+2.*  (Til,  1  ,2  I +2.*T  (  1  , 3 , 3  )  )  *  (  EN(  I  )  **2  )  * 

It  EMI  n*‘»2)+TM,2,21*EMtn*»4 

BtI,?,3)  =  tTtI,l,n-Ttl,l,21-2.*TII,3,3ll*EMtn*ENn  1*43+1111,1,21- 
lTIl,2,21+2.*Tn,3,31  1  ♦EN I  1  1*EMI  I  1**3 

- Bt l,3,31  =  tTtI ,1 ,11  +  Tt 1,2,21-2.*TI I,l,21-2.*Tir,3;3) )*tENt I  1**21* 

lIFMtn*'>21+Ttl,3,31»tEMni**4+FNtl  1**41 
BtI,2,n=BtI,l,21 
BtI,3,ll=BtI ,1,31 
Bt  I,3,21*BI 1,2,31 

C  CALCULATE  ELEMENTS_  OF  INVERSE  OF  0  MATRIX _ _ 

fo'D!  i  l ,  r,  lV*£Mt  1 1  ••2‘ 

DltI,l,21*ENin**2 
01  tI,l,31»-2.*  EMtl 1*£NI I) 

0111,2,11*0111,1,21 
0111,2,21=0111,1,11 
1)1 1  1  ,2,31=-01 1  1,1,31 
0111,3,11=0111,2,31/2. 

0111,3,21=0111,1,31/2. 

80  nl(I,3,31=EHtIl*»2-ENtIl**2 
C  CALCULATE  ELFMEMTS  OF  ALPHA  MATRIX 

AX(Il=DItI.l,l)«ALt I  1+0111,1, 21 *ATt 11 
AYt I  1=01 ( 1 ,2,1 1»AL11 1+DIt I ,2,21*ATtI  1 
90  AXYt I  1=01  I  I ,3,1 1»AL I  I  1+011 ’ ,3,2 1*ATI 1 1 
100  COMTIMOE 

UR  I TF 16,7101  I  I, 01  I  I ,1,1 1,01 11,1,21,01  I  I ,1,31 ,I ,01  I  1,2,11,01  I  I ,2,21 
1,0111,2,31,1 ,0111 ,3,11,011 1,3,2 1,01 1 1,3, 31, I  =  1,N1 
710  FtlRMATI lHl,3nx,23HFLEHFNTS  OF  01  MATRIX  /// I lOX , 1 2 , 5X , F 1 7. 8 , 5X ,E 1 
17.«,5X,F17.8  1  I 

WPITFI6,7201  t 1,B 1 1,1,1 l,Bt 1,1,2 l,BI I, 1,3 1,1, Bt 1,2,1 l,Bt 1,2,2 l,K|l  , 
12,31 ,1 ,Ht 1 ,3,1  1,B1  l,3,21,Bt 1,3,31, I  =  1,N1 
720  FORMAT! lHl,3nx,22HeLEMENTS  OF  H  MATRIX  /// I lOX , 12 , 5X ,E 17 .H , 5X , E 1 7 
l.B,5X,E17.fl}  5 

HR  ITF  16,74011  I,AX‘t  I  l,AYt  I1,AXYI  I1,I  =  1,N1 
740  FnRMATI IHl ,3nx,26HELEMENTS  OF  ALPHA  MATRIX  / // { lOX, 12 , 5X , b 17 . 8, 5X 
1,E17.8,5X,E17.R1  1 

C  CALCULATF  PROPERTIES  OF  N-PLY  LAMINATE  _ 

(TCalcolatf  elfhents  OF  b-bar  matrix 

BBI1,11=0. 

6Rtl,21=0. 

RBI  1 ,31=0. 
flBt2, 1  1=0. 

BHt2,21=0. _ _ _ 

- BB(2,3)=0.  ■■ 

BBI3,1 1=0. 

BBI3,21=0. 

BBI3,31=0. 

120  00  700  I=1,M 

OBl=Ht 1,1, 1 1*HI I  1 _ 

b8tl,ll-BHll,ll+W<T 

DB2>BII,l,21*HtIl 


o  o 


()fl3=H(  I  ) 

efi(  1,3)»BB( l,3l+0fl3 
140  0B4=H( I ,2f 1 )*H( I ) 

BR(2,  1  )=BB(2,  U+0B4 
DB5=R( 1 ,2,21*H(I  1 
Bfi(2,2)=B’.(2f2)+DB5 
nB6=B(  1 'HI n 
RB(2,3)=fil'(2,3)+0R6 
DB7rH( I ,3, 1 ) *H( ;  ) 

160  BB(3,1 )=BH(3,l)+nR7 
I  ,3,2)  >H(n 
BB(  ^,?)=hB(3,2)+llBH 
IIBq  =  B(  I  ,  3f  3)<-H(  I  ) 

BR(3,3)=BH(3,3)+0Bq 
200  CnMlINUh 

WR  1  TK6,  760)  RB(  1,  1  )  ,HH(  1,2  l,BB(l  ,3)  ,UB(?t  1  )  ,  BB  ( 2 , 2  ) ,  BH  (  2 , 3  ) ,  BB  (  3 , 
U),BB(3,2),HB(3,3) 

“  76  0  FORMAT! lHl,30X,26HFLeMFNTS  OF  b  BAR  MATRIX  // / (  1 7X , F  1 7 . 8 , 6X , h 1 7 . 6 
1,5X,F17.P1  ) 

C  CALCOLATF  FLFMF'iTS  UF  C  MATRIX 
CX  =  0. 

no  300  I=1,N 

0r.X  =  (AX(n*B(l,l,l)-*-AY(I)*B(l,l,2)  +AXY  (l)'=B(l,l,3))*Ht!) 

CX=f,X+»CX 
300  CONTIMlie 
CY=0. 

on  400  I=1,M 

OCY=(  AX(  I  )*B(  1,2,1  )+AY(  1  )*B(  I,2,2)  +  AXYI  n#BI  1 ,2,3)  )-»HI  I  ) 

CY=CY+nCY 
400  CONTINOF 
CXY=n. 

on  soo  l=l,M 

OCXY  =  (  AX  (  !  )*B(  1,3, 1 ) -FAY!  I  )»B(  1,3,2  )+AXY(  H'-B!  1 , 3,3)  )'-'H(  I  :• 

CXYiT.XY  +  n(,XY 

500  cf'MTiMue 

C  CALCOLATF  6LFMFMTS  UF  (OVe«St  OF  B-bAR  MATRIX 

OBBsBB!  1,1  )!>(BB(2,2)»BH(3,3)-BH(?,  3)*BB(  3,2  )  )+BB(  l,3)»(BB(2,ll‘>BB( 
13,2)-BB(2,2  )oBB(  3,1  )  ) +BI1  (  1 ,2  )=>  (  BB  !3 , 1  )  vBlM  2 , 3  ) -BiJ  (  3 , 3  )  »HB  (  2,  1  )  ) 
Bid ,1 )=(BB(2,2)*BB( 3,3)-BBI2,3)«BB(3,2) )/OBF 
520  HI  (  1 ,2)=-(BB(  l,2)*HB(3,3)-BHd  ,3)*BBI3,?)  )/l)HB 

Bl(l,3)  =  (HH(l,7)*nB(2,3)-BBn,3)»BB(?,2))/l-BB  . 

HI (2, 1  )=-(BR(2, 1  )*HB(3,3)-Re(?,3)*FB(3, 1 ) )/l)bB 
HI  (2,2)  =  (I<H(  1, 1  )ABH(3,3)-BH(  I  ,3)*BH|3, 1  )  )  /f)BB 
540  BH2,3)=-(BH!  1 ,  1  )  *B B  ( 2, 3  )-HB 1 1 , 3 )  *BB  ( 2 , 1  )  ) /OBB 
BI (3,1 )=(HB(2, 1 )*BH( 3,2)-HK(2,2)*BB(3, 1 ) I/OBH 
HI  (3, 2)=-(BB{  1,1  )*I)I)(3,2)-BH(1 ,2)'BB(  3, 1  )  )/Ol)H 
560  BI(3,3)=(BB{l,l)aHB(2,2)-HB{l,2)«BIM2.l))/nBU 
C  CALCOLATF  FLFMFNTS  OF  ALPHA-BAR  MATRIX 

AXB=8I(l,l)*CX+BI(l,2)*CY+BI(l,3)*r.XY 
AYP»Bn2,  1)*CX+B  1(2,2  )*CY+B  1 12,3)  =»CXY 
580  AXYB»Bl(3,l)*CXfBI(3,2I*CV-»BII3,3I*Cxy 
HRITF(6,78())AXB,AYB,AXVH 

780  FORMAT! lHl,30X,2HHfiLeMeNTS  OF  ALPHA  BAR  MATRI X//22Xt IIHALFHA  BAR  X 
l,UX,llHALPHA  BAR  Y , U X , 12HALPHA  BAR  XY// ( 17X , 6 17.B ,5X, E 17.8 , 5X ,E I 

27. HI)  _  _ _ _ _ _ _ 

CALCOLATF  STRAINS  OEVELOPEO  IN  LAMINATE  DUE  TO  UNIFORM  TEMPERATURE 
INCREASE  l)T 
EPX=nT*AXH 
EPY«0T*AYH 
EPXY=DT'»AXYe 
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C_CALCULATF  STRSSSFS  IN  EACH  PLV _ I _ 

700  nn  Aoo  i=i,n 

SX(  t  )=Hn  ill  1)«(FPX-0T»AX(  tn'FB(ltlt2)*(EPV-0T*AVUM«^b(  1  tlt3i* 

1  (FPXY-m*AXV(  I  )  ) 

S-  I  )=H(  1  i2;i  )*(FPX-nT*AX(  m*flnt?,2)*(EPY-0T*AYn  )  )^B(  I,2f3I* 

1 (FPXY-UT»AXY(  I )  ) 

_ SXY(  n=H(  I  0,1  )«(FPX-OTsAxn  )  )*B  f  I  ,3^2  )  *J  E_PV-f)L*  AV  (  1  )  l-fHl  1 .3.3>* 

1  (FPXV-I)T*AXY(  I  )  ) 

A  (  1 1 1  ,n  =  TI  I  ,1  ,1  )»n(  1 1  l,n>T<  tfl«2  >*0(  1 12  1 1  >«T(  I  ,1|3)«0(  I  t3tl) 

A(  1 , 1 ,2)  =  7(  ti  li  I  )«()(  III  ,2)«T(  I,l,2l*OUi2f2)'»T(  1,1,3)»0(  1 13>2) 
A(I,1,3)  =  1(I  ilillPOI  I,1.3)4^T(I,1.2»*0<Ii2i3KT(!.1.3)«D(1.3i3) 
A(I,2,n  =  T(  I.2,l)«n(I,l«U-fT(l,2«2)*D(Ii2il)«Tll,2i3)«0ni3il) 
AMi2t2)  =  T(I,2,l)ao(  III, 2)*T(  1,2,2  )*Dn  1 2i2)  +  T  11,2, 3>*D(  I,3f2) 

A(  r,2,3)=7(  I,2,l)<-n(  I,l,3)VTli;2,2)*Dfl,273')^f(I,2r3)«DI  1,3,3) 
A(I,3,l)  =  T(l,3,l)-0( Iil,l)«T(I,3,2)*OII,2,l)'»T(I,3,3)»OI 1,3,1 1 
A(I,3,2)=T(Ii3,l)»l)IIil,2)+T(l,3,7)*Dll,2,2)+T(I,3,3)»0(I,3,2) 
4(1,3,31  =  1(1  ,3,1  )«()(  1  , 1,3) +T(  1,3,2  )*0(  1,2,3  )'FT  (  I  ,  3, 3  )«0('I  ,  3,3 ) 

.SL(  n=A(  1,1  ,1  )*EPX+A(  1 ,  1,2)»EPV+A{I  ,1,3)*EPXV-(T(  1 , 1 , 1  > *AL  ( I  l  +  TI  1  , 

11.2) »AT( 1) )»nT 

SKI  )=A(  1,2,1  I  *EPX+ A  (  f,2, 2  )■ *=EPV'*A  ( 1,2, 3  )*ePXY-(f(  1,2,1)  PALI  I  )*T(1, 

12.2) »AT(  I  )  )«liT 

SLT( I )=A( 1,3, 1)*EPX+A( 1,3,2 )«FPVPA) 1,3,3 )*ePXY-)T( 1,3,1 )*AL( I )♦!( 1 
l,3,?)vAT(  n)»nT 
600  CnNTIN)/l^ 

_  WKlTF(6,flOOI(I,SX(I),SY{l),SXV(l),lsl,N)  _  _ 

"aoo  F(mMAT(lHl,30X,16KTHFRMAL  S  TRF  SSF  S/'/8X  ,  7HPL  Y  VTl).  ,  8X“,7HALPHA  X,lb 
1X,7HA1,PHA  Y,  I'lXihHALPHA  XY// (  lOX  ,  1  2,5X  ,  E 1 7. 0,  $X  ,E  1 7.  B,  6X  ,  E  17. 8  ) ) 
WKlTK6,9(tO)(I,SL(I  ),ST(I),S1.T(1),1  =  1,N) 

900  FORMAT! IHI ,30X,29HTHEKMAL  STRESSES  NATURAL  AXES//BX,7HPLY  NO,,  8X, 
17HA(PHA  L,15X,7HALPHA  T , ISX , 8HALPHA  LT// ( lOX , 1 2 ,5X , £1 7. 8 , 5X , E17 .8 , 

_  25X,F17.H)) _  _ 

f.O  TO  S 
END 
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EXAMPLE  PROBLEM 


To  ili'istrate  the  use  of  AC-40,  consider  a  [0/±45/0]s  laminate  at  room 
teiqierature  with  a  AT  =  10°.  For  this  case,  on  card  A,  there  are  four  plies 
per  set,  so  N  =  4.  M  is  set  to  1  for  room  temperatures  properties,  and  DT  is 
set  to  10.0  since  AT  =  10°.  For  card  B,  only  the  angles  0.0,  +45.0,  -45.0,  and 
0.0  are  input.  The  input  data  sheets  and  resulting  output  for  this  case  are 
sho\m  at  the  end  of  this  section. 

Hie  results  from  the  output  printout  are: 
a  =  2.363  xl0"6  in./in./°F 
5y  =  7.752  xlO'6  in./in./°F 


Ply 

o 

o 

<r ^10°  F 

T j^l0°  F 

a^/10°  F 

O 

0 

Ti^.j,/10 

1 

-3.7 

-79.1 

0 

-3.7 

-79.1 

0 

2 

3.7 

79.1 

133.6 

789.8 

-137.1 

-18.9 

3 

3.7 

79.1 

-133.6 

789.8 

-137.1 

18.9 

4 

-3.7 

-79.1 

0 

-l.l 

-79.1 

0 
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NUMBER  IDENTIFICATION  DESCRIPTION  DO  NOT  KEY  PUNCH 


CLEMENTS  or  ALEHA  NATftlX 


1 

0.231999918-05 

0.106700008-04 

0.0 

2 

0.649500098-05 

0-64950009F-05 

-0. 417500008-05 

3 

0.649500098-05 

0.649500098-05 

0.417500008-05 

A 

0.231999918-05 

0.106700068H>4 

0.0 

ELCNCNTS  OE  B  BAR  MATRIX 


0.4075R5138 

06 

0.A67597508 

05 

0.0 

0. 667597508 

05 

0.123674448 

06 

0.0 

D.O 

0.0 

0.H943493R8  05 

ELEMENTS 

OF  01  MATRIX 

1 

0.100000008 

01 

0.0 

0.0 

1 

0.0 

0.100000008. 

jai 

0.0 

1  • 

0.0 

0.0 

O.lOOOUUOOfc 

01 

2 

O.SOOOOOIBE 

00 

0.500000188 

00 

-0.999999946 

00 

2 

0.500000188 

00 

0.500000188 

00 

0.999999948 

00 

2 

0.499999948 

00 

-0^499999948 

00 

0.0 

3 

0.500000188 

00 

0.500000188 

00 

0.999999948 

00 

3 

0.500000188 

00 

0.500000188 

00 

-0.999999948 

00 

3 

-0.499999946 

00 

0.499999948 

00 

0.0 

4 

0.100000008 

01 

0.0 

0.0 

4 

0.0 

0.100000008 

01 

0.0 

4 

0.0 

0.0 

0.  lOOOOOOOl- 

01 

ELEMENTS 

OF  B  MATRIX 

1 

0.300199848 

08 

0.571383568 

06 

0.0 

1 

0.571383568 

06 

0.272087508 

07 

0.0 

i 

0.0 

0.0 

0.7J0UUU008 

U6 

2  ^ 

0.917091108 

07 

0.777090808 

07 

0.682477708 

07  • 

2- 

0.777090808 

07 

0.917091108 

07 

0.682477706 

07 

2 

0.682477708 

07 

0.682477708 

07 

0.789952506 

07 

3 

0.917091108 

07 

0.777090808 

07 

-0.682477708 

07 

3 

0.777090808 

07 

0.917091108 

07 

-0.682477708 

07 

3 

“=0.6e24777ffF 

07 

-0.682477708^07 

0.789952508 

07 

4 

0.300199848 

08 
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ELEMENTS  OF  ALPHA  BAA  MATRIX 


ALPHA  BAR  X 


alpha  bar  Y 


ALPHA  BAR  X'Y 


0.236320ASE-05  0.7752U99E-05  0.0 


THiRMAL 

STRESSES 

NO. 

ALPHA  X 

ALPHA  V 

ALPHA  XV 

1 

-fi.  370138936 

01 

-0.7914A974E 

02 

0.0 

> 

0.37009277E 

01 

0. 79145 752E 

02 

0.13361493E 

03 

\ 

0.370092 77 E 

01 

0.79145752F 

02 

-0.13361493E 

03 

!» 

-0.37013893E 

01 

-0.79144974F 

02 

0.0 

THERMAL 

STRESSES  NATURAL  AXES 

NO. 

ALPHA  L 

ALPHA  T  ' 

ALPHA  LT 

-0.37014160E  01 

-0.791448066  02 

0.0 

0.78977HR1E  03 

-0,137062046  03 

-0. 18861 176E 

02 

0.789778816  03 

-0.137062046  03 

0.188611766 

02 

-0,370141606  01 

-0.791448066  02 

0,0 
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II  Supplementary  notes 


13  abstract 


This  volume  summarizes  that  portion  of  the  program  concerned  with  the  development  of 
theoretical  methods  of  analysis  for  orthotropic  laminated  composite  structure.  An 
appaidix  to  this  volume  presents  a  number  of  con5)uter  programs  that  were  generated  in 
conjunction  with  these  theoretical  developments.  The  main  bod^  of  this  volume  is 
comprised  of  independent  sections  devoted  to  the  following  topics:  extensional  and 
flexural  elastic  constants,  flat  plates,  orthotropic  sandwich  plates  and  shells 
stiffened  skin  construction,  minimum  weight  analysis,  thermoelastic  relationships, 
and  joint  analysis. 

The  appendix  presents,  for  each  computer  program,  an  abstract,  data  input  instructions, 
output  description,  and  a  typical  example  problem.  The  listed  programs  cover  the 
following  subjects:  calculation  of  laminate  extaisional  and  flexural  elastic  constants; 
honeycomb  sandwich  panels  under  in-plane  biaxial  loading,  in-plane  shear  loading,  normal 
pressure,  and  certain  ccmibined  loadings;  hat-  and  Z-stiffened  panel  stability  and 
stiffiiess;  and  laminate  thermoelastic  analysis. 
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